Mechanics of Systems of Affine Bodies. 
Geometric Foundations and Applications in 
Dynamics of Structured Media 

J. J. Slawianowski, V. Kovalchuk, A. Martens, 
B. Golubowska and E. E. Rozko 
Institute of Fundamental Technological Research, 
Polish Academy of Sciences, 
5 B , Pawihskiego str., 02-106 Warsaw, Poland 
e-mails: jslawian@ippt.gov.pl, vkoval@ippt.gov.pl, 
amartens@ippt.gov.pl, bgolub@ippt.gov.pl, erozko@ippt.gov.pl 

November 29, 2010 

Abstract 

Discussed are geometric structures underlying analytical mechanics of 
systems of affine bodies. Presented is detailed algebraic and geometric 
analysis of concepts like mutual deformation tensors and their invariants. 
Problems of affine invariance and of its interplay with the usual Euclidean 
invariance are reviewed. This analysis was motivated by mechanics of 
affine (homogeneously deformable) bodies, nevertheless, it is also relevant 
for the theory of unconstrained continua and discrete media. Postulated 
are some models where the dynamics of elastic vibrations is encoded not 
only in potential energy (sometimes even not at all) but also (sometimes 
first of all) in appropriately chosen models of kinetic energy (metric tensor 
on the configuration space), like in Maupertuis principle. Physically the 
models may be applied in structured discrete media, molecular crystals, 
fullerens, and even in description of astrophysical objects. Continuous 
limit of our affine-multibody theory is expected to provide a new class of 
micromorphic media. 



Introduction 

In our former papers, cf. e.g. [20l [21] [22] and references therein, we discussed 
in various aspects the mechanics of afhnely-rigid bodies, i.e., bodies rigid in the 
sense of affine geometry. Being interesting in themselves even from the purely 
academic point of view of analytical mechanics, they were also thought on as 
a useful model of internal and collective degrees of freedom of various realistic 
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systems. Everything done in the mentioned papers concerned a single afhnc 
body (see Fig. 1). 




Fig. 1. AfBne body, degrees of freedom. 

It is interesting, however, to investigate the mechanics of systems of such bodies. 
Certain quite realistic physical applications are obvious. For example, in various 
problems of molecular dynamics the most relevant degrees of freedom are spatial 
translations, rotations and homogeneous deformations. These are so-to-speak 
the leading modes of motion of a single molecule. In molecular crystals one 
deals with aggregates, systems of such objects (see Fig. 2). 




Fig. 2. Molecular crystal or granular medium as a system of spheres or 

spheroids. 
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Atoms in a single molecule are strongly coupled to each other, whereas the mu- 
tual coupling between separate molecules is relatively weak. So, molecules have 
a well-defined individuality and the model of a system of mutually interacting 
afhnely-rigid bodies looks physically reasonable and viable. In the long-wave 
approximation one obtains asymptotically the micromorphic continuum in the 
sense of Eringen [31 [7J H] ■ This is the continuum of infinitesimal affinely-rigid 
bodies, consisting of material points with attached linear frames representing 
internal degrees of freedom, i.e., microrotations and microdeformations if using 
the terms accepted in theory of micropolar and micromorphic media. In the 
droplet model of nuclei the collective phenomena are also sometimes described 
in terms of affine degrees of freedom [5] ; the nuclei-nuclei interaction also may 
be interpreted in terms of systems of affine grains-droplets. There are vibrat- 
ing stars and other astrophysical objects. As shown by Riemann, Dedekind, 
Chandrasekhar and others, their relevant, leading, degrees of freedom in a good 
approximation are also ruled by the affine group. So, for example, a double 
star, at least when its constituents are placed relatively near to each other, is a 
system of mutually interacting affine bodies 



In all those examples (molecular crystals are perhaps most suggestive, but amor- 
phous media also provide a good example, cf. Fig. 3) the mutual displacements 
result in stresses and the mutual rotations and deformations lead to hyper- 
stresses (see Fig. 4). There is also another application, strongly related to the 
finite element method. As known, this method was elaborated primarily for the 
computational, numerical purposes. The material body is "triangulated" like in 
geodesy and represented as an aggregate of simplexes which are so small that 
in a good approximation they may be considered as homogeneously deformable, 
i.e., as affine bodies (according to the very idea of differential calculus, every 
sufficiently regular configuration may be in sufficiently small regions approxi- 




Fig. 3. Densely packed grains: microcontinua? 
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mated by affine ones). The boundary surface of the body may be then literally 
triangulated into two-dimensional, i.e., planar affine bodies, the surface finite 
elements. In traditional applications all these simplexes were used as elements 
of the mesh for the purely numerical solving of partial differential equations of 
elasticity (or viscoelasticity, plasticity, etc.). But nowadays another approach 
became fashionable, namely one combining in an inventive way rigorous ana- 
lytical, qualitative and numerical methods [151 116j . By introducing the simplex 
mesh one represents the body as an aggregate of mutually interacting small 
affine bodies. 




Fig. 4. Grains of the regular lattice become translated, rotated and deformed. 

The continuous dynamical problem based on partial differential equations is 
then approximated by a discrete one, ruled by a finite system of ordinary dif- 
ferential equations with time as the only independent variable [12] - Obviously, 
this reduction to analytical mechanics with a finite number of degrees of free- 
dom is an approximation, the better one, the mesh is more dense and inventive. 
But one can think about combining the well-elaborated methods of the theory 
of finite-dimensional dynamical systems, both rigorously analytical and qual- 
itative ones, with the numerical techniques based on the discretization. The 
procedure becomes then more reliable and the danger of computer artefacts is 
remarkably smaller than when dealing with the purely mesh-based numerical 
treatment of dynamical partial equations [T51|T7] . There is also another problem 
on a more fundamental level. Namely, one can think that in certain situations 
the simplicial-triangulation modes of motion may be more adequate than, e.g., 
the usual phononic modes corresponding to the Fourier expansion onto plane 
waves. In any case this is some reasonable alternative approach. Incidentally, 
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when we go with the size of simplicial elements to the nano-scale, the quan- 
tum treatment becomes necessary, just like in the traditional phonon-modes 
description. 

Even if assumed models of mutual interactions (described, e.g., by binary 
potential) are formally simple, it is rather a rule than exception that, both com- 
putationally and qualitatively, the mechanics of a system of material points is 
much more complicated than the mechanics of a single material point moving 
in some fixed external field. The situation becomes even incomparatively more 
difficult when our elementary entities are no longer material points but struc- 
tured objects, e.g., rigid bodies or affinely-rigid bodies. This fact is well known 
even in relatively simple systems of mutually coupled gyroscopes. The nontriv- 
ial geometry and topology of " internal" degrees of freedom brings about some 
completely new problems and difficulties, even on the very fundamental level of 
invariance principle. Although the objects move in the flat affine or Euclidean 
space (here we consider only such a situation as realistic) , the full configuration 
space is non-Euclidean. When we deal with systems, this leads to really serious 
complications like difficulty if not impossibility at all of defining correctly the 
internal motion "as a whole" (some "average" internal motion of the system) 
and the "relative" internal motions. This is similar and as a matter of fact ge- 
ometrically strongly related to the impossibility of the correct definition of the 
centre of mass for a system of material points in a curved (Riemannian) space. 

There is also another important novelty when we pass from mechanics of a 
single object moving in a Lie group as the configuration space (Euclidean, affine, 
orthogonal, linear) to a system, i.e., when the configuration space becomes the 
Cartesian product of several copies of such a group. For many reasons, both 
computationally practical and fundamentally theoretical, particularly interest- 
ing are dynamical models invariant under the group underlying geometry of 
degrees of freedom. If one deals with a single object and its configuration space 
is just a Lie group, then the only invariant models are geodetic ones, i.e., with 
vanishing (constant) potentials; Lagrangian is then simply given by the kinetic 
energy (its metric is left- or right-invariant under regular translations). The 
reason is obviously that the only functions on group invariant under regular 
translations (it does not matter left or right) are constant ones. However, if one 
deals with the system of such objects, i.e., if the configuration space is some 
Cartesian-product-power of the underlying group, then of course there exist 
nontrivial potentials, more generally — nontrivial forces, invariant under this 
group (left or right acting). This is a very simple algebraic fact. If G is the 
mentioned group and we consider, e.g., a two-body system with the configura- 
tion space G 2 = G x G, then obviously any function F : G x G — > Y of the 
form 



is invariant under simultaneous left regular translations of its arguments: 




(1) 



{91,92) >->■ {ggi,gg 2 ) 



geG 



(2) 
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where Y is an arbitrary set and / : G — > Y is an arbitrary F-valued function 
on G. In particular, F, f may be real-valued potential functions. 
Similarly, F given by 

F (91,92) = f (52.gr 1 ) (3) 

is invariant under simultaneous right regular translations: 

(91,92)1-* (919,929), g^G. (4) 

For some special choices of / one can obtain "potentials" invariant simultane- 
ously both under left and right regular translations: 

(51,32) i-> (kgil,kg 2 l) , k,leG. (5) 

An interesting example is provided by characters of finite-dimensional rep- 
resentations of G. Let T : G — >• GL(X) be such a representation; X is some 
finite-dimensional linear space and 

T(g 1 g 2 )=T{g l )T(g 2 ), T(e) = id x , (6) 

where e G G is the neutral element (identity) of the group G. Let / : G — > M. 
be a character of T: 

f(g):=TrT(g). (7) 

It is obvious that 

/ (km) = f (mk) , F(k,m) = F(m,k) (8) 

and the function F is invariant both under @ and ^j. 

All these concepts become particularly intuitive when we consider systems of 
rigid bodies, when G — SO (3, R), the group of rotations in R 3 , or its "quantum- 
spinorial" universal covering G = SU(2). The mutual (relative) rotational de- 
flections of such "rigid bodies" as molecules in a molecular crystal, or grains 
in a structured medium, usually modify the potential energy of the system 
and this potential energy function often exhibits the symmetries as described 
above. Formally the same problem may be investigated in n dimensions, when 
G = SO (n,R) or its "quantum-spinorial" universal covering G = Spin (n). 

Incidentally, let us observe that the same properties are characteristic for 
the functions of the type: 

/( 5 ):=Tr(T(/))=Tr(T( 5 ) fe ), k e Z. (9) 

If G is non-Abelian, then, obviously, the mapping 

G 3 g^T(g k ) = T(g) k eGL(X) (10) 

is not any longer the representation of G, nevertheless, the property (|8|) holds. 
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The situation becomes more complicated when G = GL (3, K) (generally 
GL(n,R)), i.e., when we consider the system of affinely-rigid, i.e., homoge- 
neously deformable bodies. As mentioned, systems with such degrees of freedom 
may be good models of molecular crystals, granular media, and other physical 
objects. There are however some delicate problems concerning the possibility 
and physical status of dynamical affine invariance. This problem was discussed 
in many our former papers [20j [21] [22] ■ We discussed there a single affine body, 
therefore the study of affine invariance was restricted to the structure of kinetic 
energy. And even at this stage the problem of affine invariance was a delicate 
matter. Namely, according to the common belief, all non-relativistic physical 
theories contain the Euclidean spatial metric tensor g as a kind of " controlling 
parameter" . Specially-relativistic theories contain the pseudo-Euclidean metric 
of the normal-hyperbolic signature; formally its status is there exactly like in 
non-relativistic physics. The metric tensor is built up into the structure of ki- 
netic energy and other physical quantities. In this way the assumed kinematical 
affine symmetry of degrees of freedom is dynamically broken and restricted to 
Euclidean subgroup (pseudo-Euclidean in specially-relativistic theories). But 
there is at least an academically-motivated interest in following the pattern of 
mechanics of rigid bodies and their systems. This motivates the search of mod- 
els of kinetic energy invariant under the left/right action of the affine group. 
Such models exist [201 (HI 121] and are mathematically interesting. But of course 
the natural question appears as to their physical utility. In mentioned papers 
certain arguments were presented. The expected application fields usually have 
to do with situations in condensed matter, microstructure, and in all phenom- 
ena where the external influences become secondary in comparison with internal 
factors. Of course, one could try from the very beginning to think in the fol- 
lowing way: The primary concept is that of kinetic energy or more generally 
of the kinetic term of Lagrangian. Metric-like tensors are byproducts of those 
kinetic terms, e.g., as configuration-dependent coefficients of quadratic forms of 
generalized velocities. Of course, one can reproach against such a reasoning on 
the basis of General Relativity, i.e., relativistic gravitation theory. This theory 
tells us that the metric tensor is a potential of the gravitational field and its spe- 
cial (pseudo-Euclidean, i.e., fiat) form is simply something like the non-excited 
"ground state" of this field. So, one could argue, it is a primary physical factor 
and the structure of kinetic energies is its byproduct. But nevertheless one can 
admit possibility of situations where the internal factors are more essential, e.g., 
in defect theory. Moreover, one knows from the solid-state physics the situations 
where really the influence of "physical" metric on the inertial properties of bod- 
ies is at least remarkably reduced. For example, the effective inertia of electrons 
moving in crystals is characterized by the effective mass tensor which differs 
from the "physical" metric and is some byproduct of complicated interactions 
in the crystalline medium. On the opposite scale of physical phenomena, in as- 
trophysics, one can also expect situations where the configuration and internal 
motion of highly condensed objects like, e.g., neutron stars is a dominant factor 
and perhaps in a good approximation to rigorous generally-relativistic consid- 
erations (based on Einstein equations) the surrounding metric field (which in 
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turn influences the star dynamics) may be expected to be a simple algebraic 
byproduct of internal mechanical phenomena. 

When systems of mutually interacting affine bodies are studied, some mod- 
els of potential energy (usually binary ones) must be used, either somehow de- 
rived from the "micromodel" (usually in a difficult way) or phenomenologically 
postulated on the basis of some rough qualitative assumptions and symmetry 
principles. This motivates what is done below in the first section. Namely, the 
concept of relative (mutual) deformation tensors of pairs of affine bodies is dis- 
cussed in details and the scalar invariants built of such tensors are constructed. 
They are expected to be basic scalar arguments of binary potentials of mutual 
interbody interactions. There is an essential novelty in comparison to deforma- 
tion scalars of single affine bodies. Namely, there exist affincly-invariant scalars 
of mutual deformations. Because of this one can consider the hierarchy of in- 
teraction models according to their invariance groups, from Euclidean to affine 
ones. There is again an open question of applicability, nevertheless, the models 
are very interesting at least from the point of view of pure mathematics and 
rational mechanics. The analysis performed in the first section might seem too 
detailed. However, without a well-defined geometric and algebraic background, 
i.e., when relying merely on analytical matrix calculus, it is almost impossible 
to avoid misunderstandings and quite real mistakes. 

1 Linear Algebra of Mutual Deformation Ten- 
sors. Transformation Rules and Invariants 

We must begin with some algebraic and geometric formalism. 

Let U, V be finite-dimensional linear spaces and U*,V* — their duals (linear 
spaces of linear functions on U, V). We deal almost exclusively with real linear 
spaces, nevertheless, all the concepts quoted below are valid as well when U, V 
are linear spaces over the complex field C. In the complex case some additional 
concepts and structures appear (e.g., the complex-conjugate spaces, antidual 
spaces, etc.), however, we will not use them. The set of linear mappings from 
U into V will be denoted by L (U, V). Obviously, L (U, V) carries the canonical 
linear structure induced pointwisely from the target space V. The symbols U* , 
V* are evident abbreviations for L(t/,R), L(V, M.) (L(Z7,C) and L(V,C) in the 
complex case). The manifold of linear monomorphisms (injections) from U to V 
will be denoted by LM (U,V); obviously, it is an open submanifold of L (U, V) 
and is non-empty only if 

m := dim U < n := dim V. (11) 

The manifold of linear epimorphisms (projections) of U onto V will be denoted 
by LE (U, V). It is also an open submanifold of L (U, V) and is non-empty only 
when dim U > dim V. The manifold of linear isomorphisms of U onto V will be 
denoted by LI (U, V). Obviously, 

LI (17, V) = LM (U, V) = LE (U, V) (12) 
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and it is non-empty only if 

dim U = dim V. (13) 

We shall use the obvious abbreviations L(U), L(V) for L(U,U), L(V,V) and 
the usual symbols of linear groups GL(U), GL(V) for LM(U,U) = LE(U,U), 
LM(V,V) = LE(V,V). 

It is important for us that L(/7), L(V) are associative algebras under the 
usual composition rule as a product and Lie algebras under the commutator 
operation: 

[A, B] := AB — BA. (14) 

Similarly, GL(J7), GL(V) arc Lie groups and L(U), L(V) are their Lie alge- 
bras GL'(U), GL'(V) in a canonical way; the exponential mapping being un- 
derstood literally as a power series in finite-dimensional associative algebras. 
GL + (f7), GL + (V") denote the components of unity in GL(U), GL(V), i.e., the 
connected groups of linear transformations preserving any of two possible ori- 
entations in U, V. To be more in accord with popularly used terms, it is 
GL + (U), GL + (V) that are Lie groups and L(J7), L(V) are their Lie algebras. 
The complements of GL + (U), GL + (V) to GL(Z7), GL(V) arc denoted respec- 
tively by GL~(f7), GL _ (V). They consist of orientation-reversing mappings 
and are cosets of GL + (C7), GL + (V). Obviously, determinants of GL + -, GL _ - 
mappings are respectively positive and negative. Being finite-dimensional linear 
spaces (thus Abelian Lie groups in the additive sense), U and V arc endowed 
with canonical translationally-invariant Lebesgue measures (special case of the 
Haar measure); they are unique up to a constant multiplier (normalization). The 
volume-preserving subgroups, i.e., unimodular groups Um(f7), Um(V) consist 
of linear mappings the determinants of which have the absolute value equal to 1 . 
Their special linear subgroups SL({7), SL(V) consist of orientation-preserving 
mappings, thus ones with determinants equal to one, 

SL(U) = Um(t/) p| GL + (U) (15) 

and so for V. The subgroups Um(C7), Um(V), SL(U), SL(V) are "structure- 
independent" in the sense that they do not depend on additional geometric 
objects in linear spaces. The same concerns dilatational groups 

m(U) = {Mdu : A ^ 0}, Dil(V) = {Aldy : A ^ 0}, (16) 

Dil+(f7) = {Mdu : A > 0}, Dil+(V) = {Aldy : A > 0}. (17) 

Let us also quote the orientation-reversing cosets: 

SL-(U) = Vm(U)f)GL-(U), Dil" (17) = Dil(Z7) f] GL~(U) (18) 

(negative extension factors A) and so for V. All other subgroups, at least prac- 
tically useful ones, are defined as ones preserving some fixed geometric object 
or figure in an underlying linear space. 
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Below we will often deal with a pair of linear spaces U, V; their ordered 
bases (frames) will be denoted respectively by 

E=(...,E A ,...), e= (...,€*,...)■ (19) 

The corresponding contravariant vector components are denoted as u A , v l , 

u = u A E A eU, v = v l e t e V (20) 

(summation convention applied). 

The dual co- frames in U* , V* will be denoted by 

E~ 1 = (...,E A ,...), e~ x = (...,e\...), (21) 

where obviously 

E A (E B ) = (E A ,E B ) = S A B , (ej) = <e\e,) = (22) 

And, as usual, we use the lower-case convention for covariant vectors, i.e., linear 
functions: 

f = f A E A eU* 1 p = p i e i £V*. (23) 

Therefore, as usual, 

/(«) = (/, u> = f A u A > p(v) = (p, v) = p t v\ (24) 

Matrix elements of linear mappings ip £ L (U, V) are meant in the conven- 
tion: 

i>E A = e^ A , (25) 

therefore, 

V- (u A E A ) = {^ A u A ) e u *P(uy = iP\u A . (26) 
We shall often make use of canonical identifications: 

L(U,V) ~V®U\ U**~U, etc. (27) 

Very often we must deal with various transposition concepts of linear map- 
pings. This is the main motivation for this introductory review of symbols. Let 
us begin with the well-known concept of the conjugate mapping. Namely, with 
any tp e L (U, V) there is canonically associated the mapping 

V>* G L(V*,U*) ~ U* ® V (28) 

given by 

ip* p := p o tp for any p G V* . (29) 
In terms of coordinates: 

(rp) A =p i Tp i A. (30) 

Roughly speaking, the matrix of ip* is transposed to that of tp in the usual an- 
alytical sense of matrix calculus. No additional geometric concept is needed for 
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defining the conjugation star- mapping (j!?9")) , (|30p . And no wonder, the canonical 
isomorphism of linear spaces 

L(U,V) ~ V"® U*, L(V*,U*) ~ {/* <g)V (31) 

described by the star-operation (|29l) , (|3"U)) is just the special case of the natural 
isomorphism between W <S> Z, Z <g)W (W, Z being arbitrary linear spaces). 

If U , V have the same dimension and ip is an isomorphism, i.e., tp £ LI (U, V), 
then we can define the contragradient mapping -0* £ ~LI(U* , V*) given by 

i>* wr 1 = ty- 1 )* ■ (32) 

Analytically: 

ty„/), := f A ^- XA i, (33) 

where obviously 

i>- XA i^B = S A B , r A ^- 1A j = 5*,-. (34) 

For any tp £ L(U, V), the upper-star-symbol ip* is also used to denote the 
natural pull-back mapping of the covariant tensor algebra over V onto that over 
U. Similarly, one uses 0* to denote the natural extension (push-forward) off/" to 
the contravariant tensor algebra over U (mapping it into the contravariant tensor 
algebra over V). If U, V have the same dimension and ip is an isomorphism, 
then ifj*, tp* may be extended to the isomorphisms of the total tensor algebras, 
and obviously, 0* = (tp*) in this extended sense as well. 
Obviously, for any tp £ L (U, V), <p £ L(V, W), we have 

((pip)* = tp*tp*, (fi>)* = V*^*, (35) 

i.e., respectively anti-representation and representation properties. 

Below we shall often use the natural isomorphism of linear spaces L (U, V)* 
and L(V, U). Namely, <p £ L(V, U) acts on tp £ L (U, V) as a linear functional in 
the sense of the trace formula: 

(tp, tp) := Tr(^) = Trtyxp) = ^aVi- (36) 

In particular, L(J7)*, ~L(V)* are identical with L([7), L(V). 

We were dealing as yet with amorphous linear spaces U, V without any 
additional structure. From now on we assume U, V to be endowed with some 
fixed metric tensors, i.e., symmetric non-degenerate bilinear forms, 

g £ Sym(V* <g> V*) C V* ® V* , r, £ Sym (U* ® U*) C U* <g> U* . (37) 

They are analytically described by their components gy, t\ab with respect to 
frames e, E: 

r) = r)A B E A ®E B , g = g ij e i ®e j . (38) 
Their reciprocal contravariant tensors denoted by 

g- 1 £ Sym(V ®V) C V <8 V, rf 1 £ Sym (U ® U) C U <8> U (39) 
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are analytically represented as 

V = r) AB E A ®E B , g = g i ie i ®e j , (40) 

where, obviously, 

r] AC rjcB = 6 A B , g lk g k , = «jy (41) 

In mechanical applications r/, g are positively definite, however, at this stage 
the definiteness problem is not essential. 

The pair of Euclidean spaces (U,i]), (V,g) gives rise to a richer system of 
invariantly defined geometric objects. First of all, one can define new transpo- 
sition concepts. They are metric-dependent, i.e., explicitly built of the metric 
tensors rj, g. Just here they are different and more peculiar than the purely 
affine (linear) concept of the conjugate mapping ip* defined above. Let us begin 
with coordinate-free definition and later on we quote analytical formulas used 
in calculations. 

The metric tensors 77, g give rise to the "index-lowering" operations produc- 
ing covectors from vectors: 

l v EL(U,U*), l g eL(V,V*). (42) 

They are given by: 

l v u = 7](u, •)=??(•, u) e U* , l g v = g(v,-)=g{-,v)€V* (43) 

for any u £ U, v € V. More explicitly, 

(l v u,w) = T](u,w), (l g v,z) = g(v,z) (44) 

for any w £ U, z G V. 

Due to non-singularity of 77, g (positive definiteness not yet necessary), linear 
operations l v , l g are invertible. Their inverses, i.e., "index-raising" operations 
will be denoted by 

r n = l~ x €L(U;U), r^l- 1 &L(V*,V). (45) 

The (r], g)-transpose of tp E L (U, V) is defined as the linear mapping given 

by 

L(V, U) 3 </; T ^ := r v ri g - V%- ( 46 ) 

If (i],g) are obvious from the context, one writes for brevity ip T . Analytically, 
ip T is given by its matrix [ - T ' 4 i] , where 

V A i = V AB g k ^ k B = 9i^ k BV BA - (47) 

If the metrics rj, g are strictly Euclidean, i.e., positively definite, then there 
exist orthonormal frames, i.e., such ones that 

Vab =* V AB =* s ab, gi 3 =* 9 13 =* Sij- (48) 
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In such coordinates the matrices of ip* and ip T numerically coincide. It is no 
longer the case in general coordinates. And in the case of non-definite (pseudo- 
Euclidean) metrics such coordinates do not exist at all. Let us observe that 
literally speaking, ip* and ip T belong to different linear spaces, respectively 

L(V*,U*) ~ U*®V, L(V,U) ^ U®V*. (49) 

They are mutually dual in a canonical way, however, they are not canonically 
isomorphic to each other. The isomorphism may be established only on the 
basis of metrics rj, g. 

The above constructions are valid for arbitrary dimensions of U, V; they 
need not be equal. If they are equal and if ip is an isomorphism, there exists 
the inverse mapping ip~ l £ L(V, U) and the contragradient mapping 

i^^ip*- 1 =tp- 1 * EL(U*,V*). (50) 
In certain formulas we shall need the (rj, g)-transpose of ip~ x , 

(ff^'eL^^-ygt/*. (51) 

It is given by 

(^" 1 ) T = r^- 1 *^ = l g -^- l %, (52) 

where the labels 77, g at the transposition symbol arc omitted for brevity if there 
is no misunderstanding danger. The matrix of tp~ 1T is given by 

&- 1T ) i A = VABip- 1B j g ji . (53) 

The above transposes depend explicitly on the both metric tensors rj £ U* ® 
U* and g <E V* ®V* . In certain formulas we use another type of transposes, 
depending only on one of the metrics. 

So, for any 

ip G L (U, V) ~ V <8> U* (54) 

we define its g-transpose 

ng) 4>eL(V,U*)-U*®V*: (55) 
again we omit the label g when it is clear from the context. Then is given 

by 

tV> := n 9 (56) 

and has the matrix with elements 

T^P Ai = 9ijtp J A- (57) 

It depends explicitly on g but does not depend on 77. If g is positively definite 
(Euclidean) and the frame e is orthonormal, then again the matrices of ip* and 
Tip are identical. Let us observe that obviously 

ip T = r vT ip = l^TTp (58) 
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and it is only here that the metric of U appears. 

If dimensions of U, V are identical and ip G L({7, V) is an isomorphism, then 
for tp^ 1 G L(V, U) we define the transpose 

t^- 1 := ip- 1 *^ € L (*7, V*) ~ ® [/*. (59) 
Its matrix is given by 

(t^T 1 )^ - VABip- 1B *. (60) 

Obviously, 

(V'" 1 ) T = r gT ^- x . (61) 

The transpose t(?7)' ! /' _1 depends explicitly on 77, but does not dependent on g. 

It is seen that the purely analytical meaning of the transposed matrix may 
be misleading when one does not take care to the geometric status of objects 
represented by matrices. In various contexts all the mentioned transposes occur 
in mechanics of deformable bodies. Everything said above is applicable to the 
special case of automorphisms, when U = V and 77 = g. Being mixed tensors, the 
elements of L(£/), L(V) cannot be transposed in an affinely-invariant manner, in 
particular they cannot be either symmetric or antisymmetric. Their transpose is 
always an essentially metrical concept. When dealing with the general L (U, V), 
one is faced with this fact in an even more drastic way. With an obvious 
exception, of course, when U = V* or U = V. Then respectively 

L(V*,V) ~ V®V** ~ V® V, L(V,V*) ~ V* V*, (62) 

and obviously the concepts of transposition, symmetry and skew-symmetry are 
well defined without any use of metrical concepts, thus, affinely-invariant. 

Having fixed the metric tensors 77, g in U, V we can distinguish within 
LI (U, V) the subset of isometries 

0{U,r,;V,g)cll{U,V), (63) 

consisting of isomorphisms preserving the metric structure, i.e., interrelating r], 
9- 

^eO(U, V ;V,g): V = V*9, (64) 

i.e., analytically 

r\AB = gn^AV 3 b- (65) 

For general isomorphisms of U onto V the above relationship docs not hold 
and the concept of deformation tensor appears in a natural way. 

For any ip e L1(U,V) the Green deformation tensor is defined in the usual 
way: 

G[ip] = ip*g e Sym (U* ®U*), (66) 

i.e., analytically 

G[(p\ab = gij^A^B- (67) 

Let us notice that G[ip] does not involve the tensor j] at all, so strictly speak- 
ing, the term "deformation tensor" is not correctly used here, as no relationship 
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between 77 and g occurs here at all. It would be more adequate to call G[ip] 
the "metric" of U which is cp- induced from 77. But for traditional reasons we 
follow the terms used in the Eringen book [SJ H] (although, to be honest, an 
"infinity" of other names occurs also in the literature; we stick to the Eringen 
conventions). Using the above transposition concepts we can simply write 

G[<p] = tW G L (U,U*) ~ U* <g> U*\ (68) 

as mentioned, to reduce the crowd of symbols, we refrain from using the unam- 
biguous expression T(g) l P l P, because g is assumed to be known and fixed. 
For certain purposes we use also the mixed tensor 

G[(p] = tp T <p S L(U) ~ U®U\ (69) 

analytically 

G[ip] A B =r, AC G[cp}cB=r l AC 9ijV i c<P j B. (70) 
Also here we do not use the more correct symbol 

G[p] = (^ T( "' 9 V (71) 

to avoid the "multi- floor" crowd of characters. If one uses rectilinear orthonor- 
mal coordinates in U, V, then the symbol <p T <p becomes exactly the standard 
notation used in elasticity. 

Some warning: one uses also the quantity 

G[<p] eU®U (72) 

given by the 77-raising of indices, analytically 

G[f] AB = ri Ac r) BD G[ip]cD- (73) 

It must be not confused with the contravariant reciprocal, 

G^}- 1 <£U <S)U, (74) 

where 

(G[cp}- 1 ) AC GMcB = S A B . (75) 

Both of them are used and coincide only when ip is an isometry, i.e., when 
<p G O {U,r}; V,g), reducing then to t? -1 , the contravariant reciprocal of 77. And 
of course, in this case, when there is no deformation, G[tp] itself becomes 77. 
Then it is convenient to use the deformation measure which vanishes. It is the 
Lagrange deformation tensor E[ip] € U* ® U* : 

E[<p] := \ [G[p] - 77) . (76) 

Unlike G[<p], E(ip) depends explicitly on both g and 77, so in this case the 
term "deformation tensor" is correctly and literally used. Just as previously, 
the objects with 77-raised indices are used, 

E[<p] eU®U*, E[<p}£U®U; (77) 
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analytically 

E[<p] A B = ^E^cb, E[(p] AB = r) AC T] BD E[tp]cD- (78) 

Obviously, in general something like E[(p] , E^]^ 1 need not be well denned 
and even if it accidentally happens to exist, it is rather non-useful to anything. 
One can suspect that perhaps some profit might follow from using the quantity 

£ [<p] £ U ® U* ~ L(U) (79) 

given implicitly by 

G[tp]:=exp(2£[tp]). (80) 

Obviously, for "almost isometric" ip ("small deformations") the quantities £[ip], 
E[tp] asymptotically coincide (exp(x) f=a 1 + x for x ks 0). 

Let us observe that G[ip], G[ip], G[<p], E[(p], E[ip], E[ip] exist also when ip is 
not an isomorphism, in particular, when U, V have different dimensions |13j . 

For any ip £ LI (U, V) the covariant Cauchy deformation tensor C[ip] is 
defined as follows: 

C[V>] = ^7 1 = Tp- 1 *r 1 £Sym{V* ®V*). (81) 

Analytically 

CMy = ilAB^- 1A l ^' 1B J - (82) 

Just the same remark as previously: it does not depend on g, does not relate 77 
to <?, so it is not a "true" deformation tensor. It is a "metric- like" tensor in V 
which is ^-induced by r\. Using the transposition symbol we can write that 

C[ip] = t^-\ (83) 

where again for brevity we write simply T instead of T(rf). 

Again some byproducts of C\ip] are important and widely used. Let us begin 
with the contravariant inverse C^] -1 £ Sym (V <g> V), also ^-independent one, 

C^]" 1 = ^.jj- 1 , (84) 

i.e., analytically: 

C[V>r Hi = VaVbti**. (85) 

The ^-dependent byproducts of C[i/)] denoted by C[ip] £ L(V) ~ V ® V* and 
C[ip] £ V <E) V are analytically given by 

CM 1 , = g lk C[ib] k] = g ik nAB ^- XA k ^- XB h CM** = g ik g jl C[^] kl . (86) 
Obviously, 

dty] = ^ 1T ^-\ C^]" 1 = (87) 

When using orthonormal coordinates and matrix language, the last expression 
coincides with the known analytic formula from elasticity textbooks. Just like 
the following one: 

CM-^VMO- (88) 
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Let us notice some important and delicate point. The covariant (thus, metric- 
like) Cauchy tensor C[ip] does exist only when ip is an isomorphism (thus, in 
particular, U, V have the same dimension). But the contravariant inverse is built 
algebraically of ip alone, without any reliance on ip -1 . So, strictly speaking, the 
primary quantity is just the contravariant push- forward of 77 : 

J[^]:=il>*V~\ JW=^A^ B ri AB , (89) 
i.e., using the transpose symbol again 

J[V]=V(tV0- (90) 

If U, V have the same dimension and ip is an isomorphism, then we can obtain 
the covariant Cauchy tensor as 

C[V>] = JW 1 - (91) 

For the general, not necessarily isomorphic ip, the primary quantity is J[ip] and 
it need not be reciprocal to anything (being in general singular). 

If ip is an isometry, then C[ip] = g. Just as previously, it is often convenient 
to use the deformation measure vanishing for isometries. The Euler deformation 
tensor e[ip] <E Sym(T/* ® V*) is defined as follows: 

e[^} = \{9- Cm) . (92) 

Notice the seeming difference in sign convention in comparison with E[ip]. It is 
not accidental, the both conventions are compatible with each other on the level 
of so-called infinitesimal deformation tensors. It would be, maybe, convenient 
to use the quantity 

s[ip] e L(V) ~V®V* (93) 

given implicitly by 

C[V]=exp(-2 £ [V]). (94) 

These were deformation tensors of linear mappings acting from U to V; they 
were always related to some fixed metric structures in U, V. The next important 
concept is one of the mutual deformation of two linear mappings from U to V. 

Let us begin from some digression. In affine spaces with every pair of points 
there is associated a unique translation vector relating them to each other. 
Similarly, with every ordered pair of group elements there are associated two 
ones which acting as left or right translations interrelate them. 

Let <p,ip G LI (U, V) be two arbitrary isomorphisms. Following the men- 
tioned concept of translation vector and left or right translations in group man- 
ifolds we introduce the following concepts of mutual displacements in LI (U, V): 

r[V>,<£] = V~V G GL(£/), £[V>, = yip' 1 G GL(V). (95) 

Obviously, 

YtyM =T[ip,iP}-\ Z[iP,ip) = Z[v,iP}- 1 , (96) 
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i.e., these operations are "skew-symmetric" in the group sense with respect to 
their arguments, just like the translation vector in afhne spaces. Obviously, the 
objects T, E are well define only for isomorphisms ip, ip. And they are equal to 
the group identity if and only if ip = <p. Incidentally, it may be also convenient 
to use "infinitesimal" versions of displacement objects, 7, er, where 

T[xl>M ^hhy+7, E[V>,y>] =U v +a, (97) 

vanishing when ip — tp. And again perhaps more convincing geometrically would 
be "logarithms" of T, E, i.e., a, (3 defined by 

rhM =exp («#,¥>]), EhM =exp(/3hM). (98) 

As we remember, the Green deformation tensor G[ip] and contravariant Cauchy 
tensor J[ip] do exist for the general ip, not necessarily isomorphic ones. It is 
natural to ask for their analogues for pairs of mappings ip, ip e LI (U, V). 
The two-argument analogue of the Green tensor is given by 

G[iP, p] e U* <g> U* ~ L (U,U*) , (99) 

where 

G[iP,<p]:= T i>ip (100) 
with the above sense of symbols. Analytically 

G[ip, lp]ab = gij^A^B- (101) 

Obviously, the symmetry is now replaced by 

G[iP,pf = G[p,iP], (102) 

here with the GL([/)-invariant sense of transposition (twice covariant tensor). 
Just like G[ip], G[ip,ip] depends on g but not on 77. The obvious ^-dependent 
object is given by 

G[ip,p] := ip T p e L(U) ~ U®U*, (103) 

i.e., analytically 

G[ip, <p\ A B - V AC G[iP, v]cb = v AC 9^ 1 CP j b- (104) 
Again in the metrically-transposed sense we have that 

G[<p,iP]=G[iP,p>] T . (105) 

The mutual Cauchy tensor 

C[ip,ip]eV* ®V* ~L(V,V*) (106) 
does not depend on g but is explicitly /^-dependent, 

C[iP,p] -.^Tp-'ip- 1 , (107) 
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<7[V>, <p]y = VABV~ XA ii>~ XB r (108) 



i.e., analytically 

/ '' ' ' '"Jij — 'lABf j 

As previously, the transposition of arguments is identical with the affinely- 
invariant tensor transposition 



C#,d T = Cfo,# (109) 
The mixed-tensor representation 

C[V>,^] G V (g> F* ~ L(V), (110) 

or analytically 

c[^p, = cpctyMki = g^VAB^k^j, (in) 

depends explicitly on both 77 and g, whereas the permutation of its arguments 
results in the (77, (^-metrical transposition. The inverse contravariant of C[ip, tp] 

Clip^}- 1 =ip T <p eV ®V ~L{V*,V) (112) 

is analytically given by 

C[lM _ " j =^a^bV AB - (113) 

Again we see that now this contravariant tensor is again more natural then 
the covariant one C[ip, ip], because no inverses of ip, tp are used. Therefore, we 
can use J[ip,ip], the analogue of the one-argument J[ip]. Even if ip, ip are not 
isomorphic, (even if they cannot be so because of different dimensions of U , V), 
we can define 

J[ip,ip] GV®V ~L(V*,F) (114) 

just as 

^rf^^Vs^. (115) 
This object is primary and if tp, ip are isomorphisms, then 

J[i>,<p] = C[iP,<p]-\ (116) 

Unlike this, G[ip, tp] is natural and its contravariant inverse G[ip, does exist 
only if ip, tp are isomorphisms. 

In analogy to ([69]), (|73| . (|86l) one uses also the tensors 



G[il>,<p], G\*I>M, CtyM, C[^M- (U7) 
Let us observe that if ip is an isometry, ip £ O (U, rj; V, g), then 

G[ip,Lp] = r[ip,<p}. (us) 

Similarly, if tp is an isometry, then 

C[iP,<p}=Z[ip,<p]. (119) 



19 



It is rather intuitive that in those special cases the mutual (two-argument) 
deformation tensors coincide with the group-theoretic "displacements" T, E. 

The next important problem is one concerning transformation properties of 
the above objects. The associative algebras L(/7), L(V) may be considered as 
semigroups acting on L (U, V) according to the rules: 

L(U,V)3tp h-> Atp, AeL(V), (120) 
L(U,V)3<p h-> LpB, BeL(U). (121) 

When restrict ourselves to A, B running over the groups GL(V), GL(U), 
then obviously Ll(U,V) is preserved under mappings (|120p . (| 12 1[) which form 
then the transformation groups acting transitively, freely and effectively on 
LI (U, V). Let us observe however that if dimensions of U, V are not equal, then 
GL(T^) acts also transitively on the manifold of monomorphisms LM (U, V), but 
the action of GL(U) is non-transitive [T31 [H]. One should also remember that 
the action of GL(V) x GL(U) on LM (U, V) given by 

LM (U, V) 3 ip i-> A<pB, A e GL(V), B e GL(J7) (122) 

is non-effective. For example, if dim U = dim V, the kernel of non-effectivity is 
given by the following subgroup: 

{ (Xld v , \-Hdu) : A e K\{0}} c GL(V) x GL(U). (123) 

The action of subgroups of GL(V), GL(U) on L(U,V) is also interesting in 
various problems. Obviously, the restriction of (|122[) to the subgroup SL(V) x 
SL(J7) is an effective action if both V, U are odd-dimensional. If they are even- 
dimensional, the kernel of non-effectivity is the two-element subgroup of (|123l) 
corresponding to A = ±1. If (|122j) is restricted to the group UL(V^) x UL(£7), 
where UL are unimodular subgroups of GL (determinants with the modulus 
equal one), then, obviously, the same two-element non-effectivity does exist 
independently of dimensions of U, V. 

Let us remind that the Lie algebras SL(V)', SL(U)' consist of traceless linear 
operators. 

In many applications, including elasticity and other branches of continuum 
mechanics, we are particularly interested in the restriction of (|122l) to the isom- 
etry groups products O (V,g) x O (U,rj), where obviously O (V,g) is an abbre- 
viation for 0(V, g;V, g); the same concerns 0(U,r]). These are subgroups of 
XJ(V,g), U (t/, rj), so in the case of even dimension their action through (|122|) 
has also the two-element non-effectiveness. Replacing the full orthogonal groups 
by their proper (determinant-one) subgroups 

SO(V,g) = 0(V,g)f)SL(V), SO (U^) = O (U, r?) f| SL(U), (124) 

one obtains again the effective action in (|122l) . Let us remind that the Lie alge- 
bras SO (V, g)' , SO (U,rj)' of the proper orthogonal groups consist respectively 
of the g- and ^-skew-symmetric operators, i.e., 

9 {Ax, y) = -g (x, Ay) , r\ (Bu,v) = -rj (u, Bv) , (125) 
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analytically 

A l 3 = -A, 1 := -g ik 9jl A l k , B K L = -B L K := -r) KM r, LN B N M , (126) 

Let us now discuss the transformation properties of the deformation and bi- 
deformation tensors. This will be necessary in symmetry studies of elasticity 
and other branches of continuum mechanics. 

For any linear isometry A G O (V, g) we have obviously 

G [Aip, A<p] = G[i/>, <p], G [Atp] = G[ip]; (127) 

the latter is just the special case of the first one, because obviously 

G[^]=G[# (128) 

And for an arbitrary A G GL( V) there is no concise expression for G [Aip, A(p] 
in terms of G[ip, ip] (also if ip = ip). But for any B G GL([7) the following nice 
formula holds: 

G[^B,<pB\ = B*G[i/>,(p], (!29) 

i.e., analytically 

G [ipB, ipB\ KL = G[4>, <p]mn B M K B N L . (130) 

For the Cauchy-type tensors we have so-to-speak the dual transformation rules. 
So, for any isometry B e 0(U,rj) 

C[il>B,<pB]=C[fl>,<p], i.e., J[ipB,<pB]=J[i>, V ], (131) 

and for a general B 6 GL(C7) there is no explicit expression for C [ipB, ipB] 
{J [4>B,(pB}) as an algebraic function of C[ip,ip] (^[tp,^]). But for any A G 
GL(V) the following analogue of (fT25j) holds: 

C[Aif!,Atp] = A*C[i>,ip], i.e., J[A1>,Aip]=A*J[1>M- ( 132 ) 
Analytically 

ClA^Mij =C[^c P ] kl A- lk l A- ll 3 , (133) 

i.e., 

J[A1>,M V = A\AhJ\^M U ■ (134) 

And finally, it is clear that the non- metrical mutual displacements T[ip, (f\, 
H[ip,tp] transform in the following way: 
For any A G GL(V) 

T[A>ip,Atp]=T[<ip,tp}, X[Aip,A<p] = AY,[iP,(p]A- 1 . (135) 

and for any B G GL(U) 

r[ipB,ipB] = B- 1 T[ip,ip\B, E[ipB,cpB]='E[ip,ip]. (136) 
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These are invariance and the adjoint/co-adjoint rules. 

In all physical applications, including ones in mechanics of deformable bod- 
ies, one needs scalar quantities built of tensors; their "invariants" roughly speak- 
ing. Let us quote the basic ones built of the tensors G, C, T, S. The standard 
(as a matter of fact canonical and unique) way of constructing scalars from 
mixed second-order tensors is to take the trace of their powers (as linear map- 
pings). According to the Cayley-Hamilton theorem it is sufficient to take n of 
them if the dimension of space equals n; any other will be some function of those 
ones. Obviously, there is no way to construct well-defined scalars from a sin- 
gle twice-covariant or twice-contravariant tensor (only scalar densities may be 
built of them) . What is then commonly called the invariant is a scalar quantity 
built of the pair of tensors, one of them being implicitly assumed. This fact is 
often overlooked when one works analytically in 1™ using formally the matrix 
language. 

Let us begin with the basic affine invariants 

M a [V^]=IV(r[^] Q )=Tr(E[^n, a = l,...,n, (137) 

where, obviously, X a denotes n-th algebraic power of the linear mapping X. 
These invariants are really affine in the sense that 

M a [AipB, AipB] = M a [ip, <p] , a=l,...,n, (138) 

for any A G GL(V), B G GL(Z7), i.e., they are invariant under the total group 
(|122p and do not assume any metrical or other additional structures in the linear 
spaces V, U. Let us observe that these quantities are essentially two-argument 
ones and, just like T, E themselves, they trivialize (become constant) when 
ip = if. There exist also scalars which do not degenerate in this sense, however, 
for the price that the above affine invariance is lost. They preassume the metric 
tensors g G V* ® V* , i] G U* ® U* and depend explicitly of them. They are 
given by 

KatyM :=Tr (dfatpf) , a = l,...,n. (139) 
One can easily check that 

K a [ip, <p] = Tr (C [j>, ^p) = Tr (J [if,, p) . (140) 

If tp = ip, the resulting quantities 

K«M:=£«hM (141) 
become that is called in elasticity deformation invariants. It is clear that 

K a [Ail>B,A<pB]=JCa[il>,<p], a = l,...,n, (142) 

for any A G O (V, g), B G O (U, rj), i.e., they are invariant under (|122[) restricted 
to O (V, g)xO (U, rj). Any other scalar built of (ip, ip) and showing this rotational 
invariance is a function of K, a [ip,(p], so they form a functional basis of the space 
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of rotational invariants. In mechanics of continua one uses also various other 
choices of basic invariants. 

Remark: More precisely, JC a span functionally the space of scalars invari- 
ant under O (V, j)xO (U, rj) but not under any larger subgroup of (|122[) ; let us 
observe that M. a [4>, <f] being GL(V) x GL([/)-invariant are automatically in- 
variant under also under the subgroup O (V,g) X O (U, rj), but are not functions 

Of K a [lp,<p]- 

The usual deformation invariants JC a [ip] (or any other system of n func- 
tionally independent expressions built of them) tell us, roughly speaking, how 
strongly the elastic body is strained, what are the scalar magnitudes of de- 
formation, i.e., stretchings, but tell us nothing about how is this deformation 
oriented in space (compare: velocity vector and its absolute value). The same 
concerns the two-argument quantities K, a [ip, ip] characterizing the scalar multi- 
magnitudes of mutual "displacements" between ip and ip. 

It is worth to mention about the volume orientation problems. Being finite- 
dimensional linear spaces (thus, Abelian Lie groups under the addition of vec- 
tors), U, V have canonical- up-to- normalization translationally-invariant mea- 
sures, namely the Lebesgue measures (the special case of Haar measures). This 
has nothing to do with other metric-type concepts like distances and angles. 
However, when some metrics rj G U* <£> U*, g G V* (g> V* are fixed, then there is 
a canonical normalization of those measures. The orientation problem is more 
complicated. Any linear space has two possible orientations and fixing 77, g does 
not distinguish any of them. Let us fix some orientations w, £ respectively in 
U, V. They are equivalence classes of n-forms, i.e., elements of A n U* A n V*; 
the equivalence given by multiplication by positive scalars. The manifold of 
isomorphisms LI (U, V) is non-connected. Its connected components consist- 
ing of orientation-preserving and orientation-reversing mapping will be denoted 
respectively by LI + (U, uj; V, £) and LI" (U, uj; V, £). Obviously, if U = V, these 
subsets are objectively defined, do not depend on any choice of orientation stan- 
dards and are simply given by the group GL+(V") and its coset GL~(V). The 
subgroup GL + (V) preserves any of two possible orientations in V and the coset 
GL~(V) interchanges them mutually. Obviously, LI+ (U, V), LI" (U, V) are ho- 
mogeneous spaces of GL + (V") x GL + (C7) acting through (|122[) . Also the action 
of the coset GL~ (V) x GL~ {U) preserves separately LI+ (U,V), LT (U,V), 
whereas the cosets GL + (V^) x GL"(C7), GL"(V^) x GL + (U) interchange them. 

Similarly, the isochoric groups XJL(V), SL(V), just like those for U, are 
objectively defined, i.e., do not assume any orientation or volume normalization, 
but for mappings acting between different linear spaces U, V it is no longer 
the case. Let (i, v denote the Lebesgue measures on U, V with some fixed 
normalizations. For example, they may be ones induced by the metric tensors 
T) G U* ® U*, g E V* ® V*. The symbol XJL(U, fi; V, v) will be denote the set 
of linear isomorphisms transforming the measure /i into v. Obviously, it is also 
a homogeneous space of UL(V) x UL(C/) acting through (I122p . And again we 
can distinguish the subset if isomorphisms preserving the oriented volumes: 

SL{U,fj,,u;V,u,£) := LI+ (U, u; V, £) f] UL(C7, w V, v). (143) 
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It is of course the homogeneous space of SL(V) x SL(U) acting through (|122[) . 
More precisely, it is an orbit of this transformation group acting in LI (U,V). 



2 Degrees of Freedom, Inertia and Kinematics 
of Affine Bodies 

Let us remind the basis ideas concerning afhne motion of a single body. This 
is necessary for formulating the system dynamics. We must begin with some 
formalism, just starting from repeating the notation from our previous paper 
[20l [2T1 [22] . Obviously, our objects move in the three-dimensional physical 
space; however, in many problems it is just more convenient (even from the 
technical point of view) to work in an n-dimensional space and only on the 
stage of real applications put n = 3 (sometimes n = 2,1, when dealing with 
degenerate dimension situations). So, let (M, V, — ») denote, as in all mentioned 
former papers, the physical afhne space; M is the underlying point set, V is the 
linear n-dimensional space of translations in M, and the arrow — > denotes the 
operation which to every pair of points x,y S M assigns the only translation 
vector xy £ V "originating" at x S M and "terminating" at y G M, i.e., 
carrying over x into y. All standard axioms of affine geometry are assumed 
[IO ll20ll2Tl l22]: let us remind them: 

• x% + yt + zi = 0, 

• for any y 6 M the mapping MBxi-^yxGVisa bijection of M onto V. 
These axioms imply in a standard way that 

xij + y±' = 0, xi = (144) 

for any x, y £ M. The linear space V as an additive Abelian group is canonically 
isomorphic with the translation group T(M). This group consists of mappings 

t v : M ->■ M, veV (145) 

such that for any x £ M, y = t v (x) £ M is the only point for which 

xfi = v. (146) 

This action of V on M as a transformation group is transitive, free and effective. 

As usual, we say that A : M — > M is an afhne mapping if there exists a 
linear mapping L(A) : V — > V, denoted also for obvious reasons as DA, such 
that 

A(x)A(y) = L(A)¥y: (147) 

Such mappings preserve all affine relations, i.e., they map straight lines into 
straight lines and preserve parallelism of straight lines. It is obvious that 

L(AB) = L(A)L(B), L (t v ) = ld M . (148) 
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In this way L is a homomorphism of the semigroup Af(M) of all affine map- 
pings of M into M into the semigroup L(V) of all linear mappings of V into 
V. Invertible elements of Af(M) form the group GAf(M), i.e., the group of 
amne transformations of M onto M. When restricted to GAf(M), L is a group 
homomorphism of GAf(M) onto GL(V). Translations are affine mappings; the 
only ones for which the second of (|148[) holds. T(M) is a normal subgroup of 
GAf(M). 

Fixing some reference point-origin, O E M, we can identify GAf(M) with 
the semidirect product GL(V) x s V; any mapping <f> E GAf(M) is identified 
with the pair 

. (149) 

The composition of ^-mappings is represented by the semidirect rule multipli- 
cation of indicated pairs. If some fixed affine coordinates are used, then the 
point x with coordinates x l is transformed into the point y with coordinates y l , 
where 

y i = A i j x j +v i (150) 

("linear-nonhomogeneous function" in the school language, simply the first- 
degree polynomial in affine coordinates). If cj> E GAf(M), then obviously 
det[A l j] 7^ 0; this need not be satisfied for a general <fr £ Af(M). By anal- 
ogy to linear transformations one uses the concepts of unimodular and special 
affine subgroups UAf(M), SAf(Af). They consist of such <j) E GAf(M) that 
respectively 

L[4>] £ UL(V), L[4>] £ SL(V). (151) 
Similarly (f> £ GAf + (Af), £ GAf~(M), when, respectively, 

L[4>] £ GL + (F), E GL-(F). (152) 

These subgroups or cosets are defined objectively in a bare (M, F, — >) with- 
out any additional structure like the volume or orientation standards in M 
(equivalently in V). If (M,V,—±) is structurally enriched by the metric ten- 
sor g E Sym (V* €5 V*), i.e., if we deal with the Euclidean space (M, V, — >, g), 
then the isometry subgroups OAf(M, g), SOAf(M, g) may be introduced as well 
as the coset OAf~(M, g). Their elements <j) satisfy, respectively, the following- 
conditions: 

L[4,]eO(V,g), L[4>]eSO(V,g), L[<t>] £ O" (V,g) . (153) 

The material space N is also endowed with an affine structure (JV, U, — >) 
with the same meaning of symbols as above. In "realistic" problems some 
metric tensor 77 £ U* ® U* is fixed and one deals with the Euclidean space, 
i.e., the quadruple (N, U, — >, 77). Usually we were dealing with non-degenerate 
situations when M, N had the same dimension n (U, V as linear spaces had the 
same dimension). Nevertheless, there are also models where dim TV < dimM 

EMM- 
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Affine mappings acting between different spaces are defined just like ones 
acting within the same space, i.e., preserving all affine relations. So, ip : TV — > M 
is afline if there exist a linear mapping L[ip] g L (U, V) such that 

f(x)(p(y) = L[tp]x^ (154) 

for any x, y E TV. Obviously, for the triple of afline spaces the following chain 
rule holds: 

L [ipi o tp 2 ] = L[<pi]L[tp2]. (155) 

The set of all affine mappings from TV to M will be denoted by Af (TV, M) or 
more detailly Af (TV, [/,—>; M, V, — >). It is itself an afline space in a canonical 
way; its translation space is Af (TV, M) . In the last expression the linear space 
V is interpreted as an affine space in the canonical way, namely: 

ui := v — u. (156) 

The Lagrange and Euler affine coordinates will be denoted by a K , y l respec- 
tively. They are fixed when some origins O E TV, O g M and linear frames 
(. . . , Ek, • • •), (• ■ • , ei, . . .) are distinguished, 

Oy = y t e l , oA = a K E K . (157) 

Affine transformation ip g Af (N, M) maps A g N with coordinates a K onto 
y g M with coordinates y l such that 

y i =x i + ip i K a K . (158) 

These are natural coordinates (x l , ip 1 k) on Af (N, M) induced by the choice 
of affine frames (O; . . . , e*, . . .) and (O; . . . , Ek, ■ ■ •) respectively in M and TV. 
The submanifold of affine monomorphisms of TV into M will be denoted by 
AfM(TV, M). If, as it is usually the case in our models, dim TV = dimTVf = n, 
then AfM (TV, M) becomes the manifold of affine isomorphisms Afl (TV, M). The 
L-operations map them respectively onto manifolds LM (U, V), LI (U, V). And 
in analogy to submanifolds LI + (U,co; V,t), LI~ {U,uj; V,£), subgroups GL + (U), 
GL + (V), cosets GL~(U), GL~(T^), one uses their L-co-images, i.e., submani- 
folds Afl+ (N,uj;M,£), Afl" (N,u;MJ), subgroups GAf+(TV), GAf+(M), and 
cosets GAf"(TV), GAf"(M). Similarly, the symbols SAf (TV, fj,, uj; M, v,l) and 
UAf (TV, [i; M, v) are L-co-images of SL (U, \x, uj; V, v, I) and SL (U, \i; V, v). And 
the groups SAf(TV), SAf(M), UAf(TV), UAf(M) are L-co-images of linear sub- 
groups SL(U), SL(V), UL(U), UL(y). In other words, one speaks about affine 
mappings preserving orientations, volume or both of them simultaneously. 

If metric tensors r/ e U* <S> U* , g E V* <S> V* are taken into account, i.e., 
if we distinguish orthogonal subgroups 0(U,i]), 0(V,g), SO (U, rf), SO (V, g), 
the orientation-reversing cosets 0~ (U,r)), 0~ (V, g), then in the affine groups 
GAf(TV), GAf(M) automatically the corresponding groups or cosets of Eu- 
clidean motions are distinguished just again as the L-co-images of the men- 
tioned linear subgroups and cosets. Roughly speaking, they are isomorphic 
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with the semidirect products of linear groups and translation groups in N, M. 
We shall use the symbols O (U, rj; V, g), E (N, rj; M, g) respectively for the man- 
ifold of linear isometries of {U,rj) onto (V,g) (i.e., such mappings <p £ L(U,V) 
that T) = <p*g) and for the manifold of affine isometries of (N, U,—¥, rj) onto 
(M, V, — >, g); obviously, E (N, 77; M, g) is the L-co-image of 0(U, tj; V,g). Just 
like the total AK(N, M), its submanifold E (N,rj; M, g) is non-connected and 
consists of two connected components. When orientations are fixed in N, M (in 
U, V), those components consist respectively of the orientation-preserving and 
orientation-reversing mappings (the corresponding <p = L [tp] are orientation- 
preserving or orientation-reversing linear mappings). 

The configuration space Q of the afhnely-rigid body is simply given by 
the manifold Afl (N, M) or rather by one of its two connected components. 
When orientations u, £ are fixed in N, M (in U, V), Q is identified with 
Afl + (N,lj; M,£), the submanifold of orientation-preserving mappings. Strictly 
speaking, the restriction to connected submanifolds is essential for classical con- 
tinuous bodies, not necessarily for discrete ones. And the more so in quantized 
models the non-connected configuration spaces seem to be admissible. If some 
reference point O G N is fixed in the material space, then Q may be naturally 
identified with the Cartesian product 

Q - Qtr x Qint = M x LI (17, V) ; (159) 

the labels "tr", "int" refer respectively to the translational and internal/relative 
motions. Just as previously, in traditional problems one restricts to 

Q~MxU + (U,u;V,£). (160) 

Usually (not always) the reference point O is chosen as the Lagrangian centre 
of mass. If the mass distribution in N is given by the positive, time-independent 
measure fi on N, then, by definition of the centre, 

f ofdfi(X) = 0. (161) 

And the affine Lagrange coordinates a K in N are chosen in such a way that 
a K (O) = 0, where K = l,...,n. 

Obviously, the identification (|159|) is meant in the sense that the affine map- 
ping <p : N — M corresponds to the pair (m, tp) where 

m = (j) ^0) = 0^ <p = L[4>\. (162) 

If affine coordinates a , y l are fixed in N, M then Lagrange and Euler 
coordinates are interrelated by (|I58[k x l are coordinates of the centre of mass 
00 in M, ip % K are internal generalized coordinates. This is the most natural 
choice of generalized coordinates in Q. It happens very often that the Lagrange 
coordinates a K are fixed once for all, often due to some physical conditions. 
Then the material space N is simply identified with W 1 . The manifold of linear 
isomorphisms LI (U, V) becomes identified with the manifold of linear frames in 
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V, denoted by F(V). Namely, (p e LI (U, V) is identified with e = (. . . , ex, ■ ■ ■), 
where ex = <fi£K and (...,£#-,...) is the natural basis of K™. If there is no 
danger of misunderstanding, are simply denoted by ipx- 

If affine motion is subject to additional constraints, the configuration space 
becomes some submanifold of Q. Usually we mean constraints imposed on 
internal motion. For example, if there is no deformation and the body is rigid 
in the usual metrical sense, then Qi nt becomes reduced to the manifold of linear 
isometries O (U, r\\ V,g), or rather to its connected component + (U,co; V,£). 

If the body is incompressible and orientation-preserving, Q mt is reduced to 
SLI (U, //, w; M, v, £), i.e., Q is reduced to SAfl (U, fi, w; M, v, I). 

In usual models there are two inertial objects, namely the total mass m and 
the second-order (quadrupole) moment of Lagrangian mass distribution: 



m 



J dfi(a), J AB = J a A a B d f i(a) = J BA . (163) 



The scalar m and the Lagrangian tensor J € U <S)U are constant and char- 
acterize respectively the inertia of translational and internal (relative) motions. 
Obviously, one can calculate all possible moments of this kind, 

J(k) Al '-' Ak = J a Al ■■■a AK dn{a), (164) 

however, J(l) = according to our convention concerning coordinates a K (their 
origin in U is placed at the Lagrangian centre of mass distribution /x) and high- 
order moments J(k), k > 2, are irrelevant for the affine motion. The second- 
order moment is algebraically equivalent to the one that is usually called the 
co-moving inertial tensor. 

Being explicitly dependent on the configuration (f>, thus, on time, the Euler 
moments, i.e., multipolcs of the transported measure [i^ with respect to the 
spatial coordinates y % are non-useful as inertial characteristics. This fact is 
known from the (metrically) rigid body mechanics. 

As usual, the motion is described by the time-dependence of the configura- 
tion (j), i.e., by curves 

WL3t^ 4>(t) e Q, (165) 
or equivalcntly, in terms of generalized coordinates 

I3(4i'((), M.3t^(p l K (t), R3i4eV(t). (166) 

Generalized velocities are denoted by v l , V 1 k- 

■■«>-£■ ^o-^. <-) 

Obviously, v £ V, V € L (U, V). Often one uses also Lie-algebraic quantities 
n e L(V) ~ GL(V)', n e L(U) ~ GL(U)': 

n = ^<p-\ n = L P - 1 ^f = ip- 1 n l p, (168) 
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i.e., analytically 

n i dip 1 A ia PjA ,„-lA d<f l B fiaQ-\ 

" 3 ■• = —j-j-P 3 ' il B = f ^ 169 ^ 

Roughly speaking, f2 is the co-moving representation of f2. If the motion is met- 
rically rigid, the above objects become respectively g- and 77-skew-symmetric, 
i.e., elements of Lie algebras SO (V, g)' , SO (U, r])' of orthogonal groups SO (V, g), 
SO (I/,??): 

n'j = -fi/ = -g^n^, o A B - - -vbkv al ^ k l- (170) 

This is then the angular velocity respectively in laboratory and co-moving rep- 
resentations. For the general affine motion one can introduce the g- and 77-skew- 
symmetric parts of Cl, O: 

w'j := \ (rtj - n/) , ^ B := \ (n A B n B A ) . (in) 

This notation is not very happy, because if <p is not an isometry, then w A b are 
not co-moving components of oj: 

u A B ± ip-^jw'jtfB. (172) 
Obviously, if motion is metrically rigid, ip G LI (U, 77; V, (7), then 

u i j = n i j , w A B = & A B = v~ 1A iU i y i B. (173) 

Similarly, one introduces distortional velocities as symmetric parts of fi-objects: 

d*j := i (n*,- + n/) , d A B : = 1 (n A s + h B A ) . (174) 

In rigid motion they vanish, of course, and in deformative motion, d A b again 
are not co-moving components of d l j. Let us introduce d e U* <S> U* given by 

d A B=VAcd c B - (175) 

One can show that 

-7 IdGAB /T7G\ 

iAB = 2-*r- (176) 

The quantity occurs in Eringen's papers under the name "gyration" . It is 
directly related to the gradient of the Euler velocity field. Indeed, if our affincly- 
dcformable continuum fills the whole physical space M, then the velocity field 
in M is given by 

V (yY = v* + fi*,- (y> - x>) = —7- h STj & - x>) . (177) 
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In certain formulas it is convenient to use the co-moving representation of 
translational velocity, v 6 U : 

d=(p- 1 v, v A =ip- 1A i v i , (178) 

Transformations (|122l) affect the affine velocities fi, f2 according to the rule: 

tp h- > AtpB : fl^AflA- 1 , h^B^QB (179) 

for any A e GL(V), B € GL(U), and ip E LI (J7, V). 

Similarly, the translational velocity transforms then as follows: 

ip h-> A<^B : w i— > Aw, u i y B~ x v. (180) 

Let us mention some important point. The holonomic rigid body constraints 
may be alternatively described in an apparently non-holonomic form, just by 
stating that affine velocity is skew-symmetric in the g-sense: 

n l j + n/ = n l 3 + g ia g jb n b a = o (i8i) 

Expression on the left-hand side is a combination of generalized velocities V % a = 
d(p l A/dt with (^-dependent coefficients. Nevertheless, the corresponding Pfaff 
problem is integrable and the family of integral surfaces is labelled by metrics 
i] E U* <E) U* and consists of submanifolds O (U, r\\ V, g) C LI (U, V). Obviously, 
in realistic problems r\ is fixed once for all and only one of all possible integral 
"leaves" is used as gyroscopic constraints. 

Alternatively, one can use the co-moving nonholonomic description in which 
Q is 77-skew-symmetric: 

n A B + h B A = n A B + v Ac r lBD n D c = Q (182) 

Again the Pfaff problem is completely integrable and the corresponding inte- 
gral foliation is labelled by metrics g £ V* <E) V* and consists again of sub- 
manifolds O (U, 77; V, g) C Ll(U,V). In "true" physics g is fixed once for all, 
pre-established, and one obtains the same results as that based on (|181[) . 

In certain traditional books, constraints described according to the pattern 
like (|181[) . (|182p are called " semi- holonomic" . Strictly speaking, they give some 
integral foliation and no transition between different leaves is possible. Some 
particular leaf is fixed by physical conditions and this leaf just represents true 
holonomic constraints. 

In a quite similar way the incompressibility constraints may be described in 
the following semi-holonomic form: 

Trf7 = Trfi = (183) 

The integral foliation is labelled by the possible volume standards fi, v in 
N, M and consists of submanifolds UL (U, fi; V, v) of volume-preserving isomor- 
phisms. When fi, v are fixed (there is a canonical choice in Euclidean spaces 
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(U,rf), (V,g)), then some particular leaf is fixed as holonomic constraints of 
incompressible motions. 

In the above examples constraints (|181[) . (I182[) were semi- holonomic because 
the subspaces of matrices (|181l) . (|182l) . (|183l) were commutator Lie algebras. 

Let us remind that the subspaces of g- and 77-symmetric linear mappings, 
i.e., ones satisfying 

Cl^ - Cl/ = - g la g ]b Cl b a = 0, (184) 

h A B -n B A = Q A B -7j Ac r, BD Q D c = 0, (185) 

are not Lie subalgebras of L(V), L(J7). Moreover, they are so-to speak as 
far non-subalgebras as possible; the commutators of (g, ^-symmetric mappings 
are just (g, 77)-skew-symmetric. Because of this, constraints (I184[) , (|185[) are 
non-holonomic, again let us say, as non-holonomic as possible. Physically we 
would say that they are constraints of rotation-less motion. This is the only 
reasonable definition of the rotation-free behaviour. As a physical application 
we can realize, e.g., homogeneously-deformable suspension in a viscous fluid. 
The friction at the surface prevents the body to rotate, but it is not an obstacle 
against deformations. 

Let us now remind briefly the basic Hamiltonian concepts. The tangent 
bundle TQ of the configuration space Q may be identified simply with the 
Cartesian product 

V v = M x LI (U, V) x V x L (U, V) . (186) 

Its elements (x, ip; v, V) are quadruples consisting of the centre of mass position 
x £ M internal configuration ip £ L1(U,V), the translational velocity v £ V, 
and the generalized velocity of internal motion V £ L (U, V). But as well we can 
use a representation based on non-holonomic velocities, first of all the following 
ones: 

Vn := M x LI (U, V) x V x L(V), (187) 
V h := M x LI (U, V) x U x L(U), (188) 

where the elements are respectively (x, tp; v, CI) and ^x, (p; v, Clj involving, as we 

see, the laboratory (spatial) (v, fi)-description of velocities and their co-moving 

(material) L?, Clj -representation. 

The Hamiltonian phase space, i.e., cotangent bundle T*Q, is identified with 
the Cartesian product 

Pv := M x LI (U, V)xV* x L(V, U). (189) 

Its elements (x, f',p, P) comprise the configuration [x, ip) and canonical mo- 
menta p £ V* (translational one) and P £ L(V, U) ~ L (U, V)* (the internal 
one); one makes use here of the canonical isomorphism between L (U, V)* and 
L (V, U) given by trace formula, 

(P, V) = Tr (PV) = P A l V i A - (190) 
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And just like previously, one introduces the objects dual to affine velocities, 
namely affine spin in the laboratory and co-moving representations, £ 6 L(V), 
£ G L(C/). They are given by 

£ = tpP, £ = Pip = ip^Zip, (191) 

i.e., analytically 

^ J = V \ ] P A 3 , Y A B = P A l ^ B . (192) 

Let us notice that the inverse mapping ip~ x does not occur here, so they are 
defined globally on the total L (U, V). They are also well defined when dim U < 
dim^, but, of course, the similarity relationship in the last expression of (|191[) 
does not exist then, because it is based on the existence of ip^ 1 . So, in such a 
situation, £ is not a co-moving representation of £. 
The manifolds of quadruples 

(x,(p;p,P), (x,(p;p,T,) , (x,(p;p,T^ , (193) 

i.e., various models of the Hamiltonian phase space T*Q will be denoted respec- 
tivelyby P, P s , Pg. 

Affine spin is in a sense conjugate momentum of non-holonomic velocities, 
because 

Tr(Py) = Tr (£fi) = Tr (SfiJ . (194) 

It is clear that this trace expression corresponds to the canonical isomorphism 
between L (U, V)* and L(V, U), L(U)* and L(U), L(V)* and L(V). 

Let us mention also the co-moving representation of canonical translational 
momentum, 

U 3p:=p*p = poip 1 PA=Pi<p t A- (195) 

As usual, the relationship between (v, V) and (p, P), or ones between (v, ft) and 

(p, S), yv, flj and ^p, S^j , is a dynamical concept depending on the Lagrangian 

(in usual dynamical models L = T — V(x, tp) — just on T). Obviously, (p, S) 

and ^p, Y^j are respectively Hamiltonian generators, i.e., momentum mappings 

[1 of GAf (M) and GAf(JV) acting on Q = LI (U, V) through the left and right 
superpositions, 

^4of B, (196) 

where 

(A, B) e GAf (M) x GAf(N), A e GAf(M), B e GAf (A). (197) 

Therefore, Pi, pa are respectively basic generators of spatial and material trans- 
lations. P A i are basic Hamiltonian generators of local additive translations in 
LI (U, V). And S l j, 'E A b are basic Hamiltonian generators of the left and right 
multiplicative translations (| 122[) in LI (U, V), i.e., spatial and material rotations 
and homogeneous deformations. Let us remind that just like in the case of angu- 
lar momentum one speaks about "orbital", "internal" (spin) and total angular 
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momentum, one can, just should, use their counterparts in mechanics of afhnc 
bodies. The canonical afhne "orbital" momentum with respect to some fixed 
"origin" O G M is given by 

A[0] := Ool®p = Ox ®p G V ® V* ~ L(V), (198) 
i.e., analytically, simply 

Kpf^x^. (199) 
The total affine momentum with respect to O G M is given by 

J [O] := A [O] + S. (200) 

It is a system of Hamiltonian generators of centro- affine mappings, i.e., of the 
subgroup GAf (M,0) C GAf(M) preserving O G M and acting on the left 
on configurations via (|196|) . The g-skew-symmetric parts of J{0\ : A[0], £ 
are respectively the total, orbital, and internal (spin) angular momenta. They 
generate the left (spatial) action of the Euclidean subgroup E(M, g) C Af(M) 
operating via (|196[) . In the velocity representation, via the inverse Legendre 
transformations, those all quantities have to do with the dipole moments of the 
distribution of kinematical linear momentum. The 77-skew-symmetric part of £ 
was by Dyson called "vorticity" [5 

V A B := % A B - % B A = £ A B - V ac vbd^ D c- (201) 

Remark: V A b is not the co-moving representation of spin 

:= _ = ~ 9 lk 9ji^ l k, (202) 

i.e., 

S l 3 ± <p i A V A BV- 1B j , (203) 

unless the motion is metrically rigid, ip G O (U, r/; V, g). 

Transformation rules of S, S are identical with those of f2, il (|179[) . thus,, 

ip^AtpB: S h- > AEA^ 1 , T, ^ B'^B (204) 

for any A G GL(F), B G GL(£7). 

As usual, the most convenient way of deriving equations of motion for 
Lagrangian-Hamiltonian models is to use the Poisson brackets for the basic 
dynamical quantities and write down equations in the following form: 

dF . . 

~dt ={FlH}l (2 ° 5) 

where F runs over some family of functionally independent functions on the 
phase space. Usually this family includes the above quantities. The latter ones 
satisfy the Poisson constants rules discussed in [20l [21] [22] and reflecting the 
affine group structure constants. 
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Let us remind some of them: 



{EV.E*!} = 8 i l X k j -8 k j Z i l , (206) 

{A l j,A k i} = S i iA k j -6 k j A i l , (207) 

{Tj,J k i} = S i l J k j -S k j J i l , (208) 

[£ A B ,Z C D } = 5 c B t A D - 8 a d YF b , (209) 

= 0, (210) 

\p A B ,Pc) = 6 A C PB, (211) 
{J i j,Pk} = {A i j , Pk } = 5\ Pj . (212) 

If F is any function depending only on generalized coordinates x l , <p l Ai then 

{F,^} = ^ A ^-, (213) 

= ^ (214) 
Pip 

= ^ Brr . (215) 

3 Towards the Dynamics of Systems of Affine 
Bodies 

Let us now discuss some peculiarities of the system of affine bodies. Certain 
quite new problems appear there, especially, ones concerning the dynamical 
affine invariance. In mechanics of a single affine body there exist affinely- 
invariant models of the kinetic energy, i.e., affinely-invariant (under left or right 
regular translations) metric tensors AQ.(N,M). But of course there are no 
affincly invariant models of the potential energy, because the groups GAf(TV), 
GAf (M) act transitively on Afl (N, M) and the only invariant functions are con- 
stant ones. However, potentials of mutual interactions between different affine 
bodies MAY BE affinely invariant. In any case there exist affinely-invariant 
scalar functions on the configuration space of the system of affine bodies. This 
again raises the question as to the applicability of such affine scalars as argu- 
ments of the total potential energy. Scalars of this kind were constructed above; 
take e.g., fjl3T[) . and the question arises as to what extent, if any at all, they 
may be used as arguments of the potential energy, perhaps as ones in addition 
to the well-established orthogonal invariants JC a (|139|) . 

Let us now consider the system of M affine bodies. To be not too abstract let 
us think that those bodies are molecules of molecular crystals, colliding fullerens, 
inclusions or gas bubbles in fluids, or some finite elements, e.g., in solids. 

The configuration space Q (Af) of such a system is the Cartesian product, 
Af-th Cartesian power of Q ~ M x LI(U,V), the configuration space of an 
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individual affine object: 

Q (TV) =Q M = x u Q. (216) 
Depending on what is more convenient, we can use alternative representations: 

Q (AO ~ Qtr (AO x Q int (AO - M M x LI (U, vf ~ M 1 ^ x f (V)*", (217) 

with the same provisos as previously that for classical continuous systems it is 
rather connected component of the above manifold that should be taken as the 
configuration space. However, even in the case of such objects like molecules 
this restriction may be too strong. 

In this way, configuration of systems of affine bodies (let us keep in mind 
something like molecular crystals, or finite elements to be more concrete) are 
described by arrays 

(...,x(/C),...;...,^(/C),...), (218) 

where AC = 1, ... ,7V, x (AC) E M, and ip (AC) 6 LI (U, V) or, equivalently, 

(...,x (£),?(£),...) (219) 

depending on whether we order separately translational positions and inter- 
nal/relative variables or we prefer to order one by one the total configurations 
of separate objects. Similarly, we shall use the Newton state space 

TV (AO, Pn(A0, V n (M) (220) 
and the Hamiltonian phase spaces 

iV(A0, Pn(A0, PnW- ( 221 ) 
They consist respectively of the ordered arrays of following objects: 

(x(K),<p(K);v(lC),V(lC)), (x (K) , <p (AC) ; v (K) , O (JC)) , (222) 

(x{K), V >{K);v(lQ,n('Q), (x(JC),ip(JC);p(IC),P(JC)), (223) 

(x (JC) , (JC) ;p(tC) , E (AC)) , (x (JC) , p (JC) ;p(JC) , £ (AC)) , (224) 

where AC = 1 , . . . , TV and, obviously, 

v(JC), V(JC), O(AC), v(JC), O(AC), P (JC), P (JC) , S (AC) , p (JC) , S (AC) 

(225) 

are just the previously introduced objects labelled now, in addition, by the index 
AC referring to individual "particles". 
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So, for example, 
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It is rather dull to write (the more to read) about the long list of canonical 
identifications. Nevertheless, sometimes this care about formal details (e.g., to 
make a distinction between V and 17) is convenient and prevents us from doing 
mistakes. Concerning the above particular details: there is some delicate point, 
not taken here into account, but perhaps important in certain problems. It is 
not quite clear, namely, if really always all constituents should be described with 
the use of the same material space U. Perhaps for any /C-th constituent its own 
Lagrangian space U (Kf) should be used? Here we do not go into such details 
and assume all U {IC)-s to be identical with some standard U. In any case, 
the aforementioned identification of U for a single object with W 1 (physically 
R 3 ) justifies the use of a common U, just R n , for all affinc constituents of the 
considered body. There are, however, some delicate points left. Namely, analysis 
of symmetries is an important part of our study. The groups GL(y), GL(f7) 
act in an obvious way, on the left and on the right, respectively, on the manifold 
of internal degrees of freedom Q- lnt ~ LI (U, V) ~ F(V) of a single object. By 
analogy, any A G GL(y), B G Gh(U) acts on Q int (AT) ~ LI (U, V) M separately 
on all factors of the Cartesian product. But perhaps we should take the arrays 

(...,A(/C),...)eGL(y)^, (...,£? (/C) ,...)& GL(U)^, (232) 

and let them to act on Qint (A/") — LI (U, V)^ so that each A (Kf), B (Kf) act 
separately on any ip (K) G LI (U, V)l And then the idea of taking independent, 
different U (ZC)-spaces for the /C-th objects would be reasonable. May such 
transformations have some physical meaning? It is still an open question. Let 
us notice, incidentally, that taking various transformations for various /C-s would 
be a kind of discrete gauge idea! 

Let us now discuss Hamiltonian (and also non-Hamiltonian) dynamical mod- 
els of multibody affine systems and the problem of their affine invariance. The 
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simplest and most viable models are potential ones, when Lagrangians have the 
following form: 

L = T-V, (233) 

where T is the kinetic energy quadratic in generalized velocities and V depends 
only on the configuration, i.e., on generalized coordinates. Geometrically T 
is equivalent to some Riemannian structure on the configuration space. The 
corresponding Hamiltonians are obtained in a standard way via the Legendre 
transformation and have the form of phase-space functions, 

H = T+V 1 (234) 

where T is the kinetic Hamiltonian. If in some generalized coordinates on 
the configuration space T is given by 

1 dq a da? , . 

then 

r=ir a Vp/3, (236) 

where T a p, usually depending on q M , are interpreted as covariant metric coeffi- 
cients, and T a P are coefficients of its contravariant inverse, 

r a T M/3 = S a p> ( 237 ) 

The quantities p a are canonical momenta conjugate to coordinates q a (more 
precisely, to generalized velocities q a ). The Legendre transformation for La- 
grangians (|233|) reads 

"'w'^-i- (238) 

We did not take into account forces depending on velocities, e.g., magnetic 
ones. For general models we have obviously the standard rule: 

* - w im 

Even if one admits dissipative forces like the viscous friction, the first step is 
always to start from Lagrangian model with its main kinetic term T . Roughly 
speaking, this term characterizes inertial properties of the object. Later on 
one introduces forces of variational structure, i.e., derivable from some poten- 
tial V(q) or generalized potential V (q,q) depending also on velocities dq a /dt. 
Generalized potentials may describe, e.g., magnetic or gyroscopic forces. And 
the final step, after deriving Euler-Lagrange equations corresponding to the La- 
grangian (|233|) is to introduce, either phenomenologically or on the basis of 
statistical-mechanical considerations, some essentially non-Lagrangian, dissipa- 
tive forces V^; the corresponding equations of motion have the well-known form: 

D dL dL . . 

mW-W = v - (240) 
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The problem of the main inertial term of nonrelativistic mechanics, 



was discussed in its various aspects in (20J (2TJ [22] in the theory of a single afhne 
body, certain more general aspects were also discussed by Capriz O 0] . 

It is not only theoretically interesting, but very often just computationally 
effective to use Hamiltonian formalism. One performs the Legendre transforma- 
tion (|239l) . one takes its inverse (in the "usual" mechanics it is assumed to exist), 
i.e., one expresses generalized velocities g M as functions of canonical momenta 
p v and generalized coordinates q v , and then one constructs Hamiltonian H as a 
function on the phase, i.e., on the cotangent bundle T*Q over the configuration 
space Q; analytically speaking, H depends on q^, p^, and perhaps explicitly on 
time t (if L does so) . Let us remind that H is given by the classical formula: 

QT 

H(...,g»,...;...,p, : /•(....,/"....:....</"••••;• P„ = ^j, (242) 

where E is the energy function on the tangent bundle TQ locally given by 

FiT 

E = ^~ L - (243) 

As mentioned, L, E, H may also depend explicitly on time, although in the 
above formula we do not indicate this. Although the non-dissipative Euler- 
Lagrangc equations, i.e., (|240[) with T>^ — 0, are traditionally transformed into 
the canonical Hamilton form: 

dq^ = dH_ dp, = _dH_ 

dt dp^ dt dq*' 1 ' 

in many problems, both purely theoretical and computational ones, it is more 
convenient to use not p^, but some their functions F and to write equations 
of motion in the following form: 

^ = {F,H}, (245) 

where the time differentiation on the left-hand side acts on the time dependence 
of F through g M , p^, and Poisson brackets on the right-hand side are analytically 
expressed by the usual formula: 

r „„, dF dG dF dG , nln . 

{F > G} = Wop-»-W»W (246) 

The complete system of equations of motion is obtained when the function 
F in (|245[) runs over some system of 2/ functionally independent functions on 
the phase space; / is the number of degrees of freedom, i.e., dimension of the 
configuration space. Usually one chooses as F generators of some important 
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groups of canonical transformations. An important technical tool is the follow- 
ing property of Poisson brackets, additional one to its Lie-algebraic properties: 



{f(F),G} = f(F){F,G}. (247) 

Our generalized coordinates are functions x (/C)\ <p (IC) 1 a, whereas the cor- 
responding canonical momenta are denoted by p (/C)j, P (IC) A j. Similarly, with 
every of TV objects there are associated affine velocities and affine momenta, 



n(jcy J7 n(ic) A B, s(/c)%-, %(ic) A B . (248) 

Quantities corresponding to different "particles", i.e., labelled by different IC- 
indices, have mutually vanishing Poisson brackets. Obviously, the ones for a 
given fixed particle are identical with the well-known expressions for a single 
affine body [201 EH [22] . 

Assume that there is no mentioned "discrete gauging" and all elements of 
the configuration arrays 

(...,*(£),...;...,¥>(£),...) (249) 

are transformed in the sense of ([120]) . ([120]) . (fl"22j) . (fT96]) with the use of the 
same /C-independent mappings A, B. For example, translations t v : M — t M 
and internal mappings A S LI(V), B G LI(£7) act as follows: 

(. . . , x (IC) tp (IC) ,...) ^ (. . . , t v (x (£)) A v (IC) B, . . .) , 

(250) 

where v € V, A € GL(V), B € GL(U) in a given transformation are fixed and 
do not depend on IC. Geometric meaning of the phase-space quantities 

Pu pa, £*,-, S A B , A'y, (251) 

as Hamiltonian generators implies that they are additive and that it is mean- 
ingful to introduce the total objects characterizing the multiparticle situation: 

Pi = Y, P* C^) ' PA = Y,PA (IC) , (252) 

K IC 

S*j = X] ^ (£) , = ^ t A B (IC) , (253) 

jc k 

K l 3 = ^ A*j (IC) , J* j = 53 J 1 , (IC) , (254) 

/C K 

etc. The same concerns, e.g., the corresponding doubled (/-skew-symmetric parts 
oiY, 1 j, h 1 j, J 1 j, i.e., spin, orbital angular momentum and total angular momen- 
tum. All those quantities, just like (|252j) - (|254[) generate global transformations 
in the multiparticle phase space and satisfy reasonable and geometrically well- 
motivated balance laws. Those become conservation laws if the corresponding 
dynamical models are invariant under, e.g., (|250l) . 
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Unlike (|252[) (I254[) , the corresponding composition rules for kinematical 
quantities 

v l ()C), v A {K), rtj{K), n A B (IC), (255) 

and so on are non-additive, in general complicated, and, the most important 
thing, non-objective. By this we mean that they depend in an essential way on 
the dynamical model and follow from the geometric rules (|252j) - (l254[) via the 
inverse Legendre transformation. 

What concerns the models of kinetic energy of a system, the most natural 
idea is to assume the additivity as well, just like in (|252l) - (|254[) . There is no 
logical necessity here to follow (|252[) - (l254p . where this was a consequence of first 
principles. Nevertheless, it seems to be a priori a rather natural conjecture. 

So, let 

T{K) = T (x (K) ,ip(K);i(K),tp (£)) (256) 
denote the kinetic energy of the (/C)-th particle, and 

T{K) = T (x (/C) , v (£) ; p (£) , V (JC)) (257) 

be its canonical representation in the /C-th phase space, i.e., the image of T (/C) 
under the corresponding Legendre transformation. According to the additive 
model we would have 

T = ^T(/C), r=£T(/C). (258) 

K K 

The key problem is then the structure of one-particle contributions. In this 
treatment nothing was said as yet about this important problem; we merely 
referred to our models presented in earlier papers 20, 2"Tjl22], Below we remind 
them briefly. Before doing this let us mention, however, what would be the 
simplest idea competitive to (|258[) . Obviously, it would be a binary one, 

t = \y,t{^l) 1 r=^r(/c,£), (259) 

completed by the symmetry assumption: 

T(fC,C) = T(jC,K.), T{K.,C) = T{C,K). (260) 

The additive model (|258l) is the special case of (|259[) . namely one correspond- 
ing to 

T(K ) C)=T(K)8 K c ) T{K,C) = T{JC)8 KC . (261) 

The hypothetical expression (|259p in general need not be diagonal with re- 
spect to the particle labels. 

Incidentally, let us mention some important interpretation problems concern- 
ing the choice between postulates (|258[) . (|259[) . In a sense this may be perhaps 
a semantical problem. Namely, in fact it seems rather natural to understand 
the "literal" kinetic energy in the additive form (|258[) . But it may be also rea- 
sonable to admit expressions like (|259[) as terms of Lagrangians, interpreting, 
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however, the non-diagonal terms with /C ^ £ ("crossing terms", "interference 
terms") as contributions describing some kind of interaction between different 
bodies, based on the coupling of velocities. Let us mention that such terms may 
appear in certain problems of the dynamics of systems of rigid bodies (gyro- 
scopic systems). In a sense, the couplings between spins resemble this structure 
of interactions. 

Now let us remind our earlier ideas concerning the kinetic energy of a single 
affine body [2D1 HH . 

The traditional "d'Alembert" model, closest to our elementary intuitions, is 
derived from the assumption that our affine body is an aggregate of material 
points. Using intuitive physical terms we may realize, e.g., a "molecule" consist- 
ing of "atoms" . Then, following the standard procedure of analytical mechanics, 
we perform the summation of elementary kinetic energies of constituents, then 
substitute affine constraints (homogeneous deformations), and finally obtain the 
classical formula 

T = T tt + T int = ^ a,,,",-- + X - aij V l A V^ B J AB . (262) 

The meaning of symbols is like in (|158l) . (|163p . thus, m is the total mass and J 
is the co-moving, thus, constant, tensor of inertia. The one-particle generalized 
coordinates x z , ip z a refer respectively to the spatial position of the center of 
mass and to the internal/relative motion. The symbols v i , V 1 a denote the 
corresponding generalized velocities, respectively, translational and internal ones 
(|167p . The above expression (|262j) is invariant under the spatial action of the 
group O (V,g) and the material action of O (U, J^ 1 )', we mean transformations 
in the sense of (|122[) . Obviously, it is also invariant under translations in M and 
TV. Particularly interesting is the special case of inertially isotropic body, when 

J KL = Ir) KL , (263) 

where I is some positive constant. Then O (U, J~ x ) — 0(U,rj), i.e., kinetic 
energy is ^-isotropic. 

It is interesting in itself but also suggestive and inspiring in the perspective 
of certain alternative models of T to rewrite (|262[) in other equivalent forms, 
namely, 

Ttr = yfltf«V = ^G[v]abv A v b , (264) 
= aV 3 bJ AB = \G\ip\ K jSl K A ^l h bJ AB 

= ^gijWWkJ&f, (265) 

where £l K \, v K are given respectively by (|169[) . (|178l) . and J[<p] is the 
spatial, thus, (^-dependent, inertial tensor: 

J[ip] kl =^ k A ip l B J AB . (266) 
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Let us remind that J\f] may be expressed in terms of the configuration 
<j) : N -¥ M as follows: 

J[<P\ kl = J (y k - * k ) (y l - x l ) d H ( y ) t (267) 

where fi^ is the ^-transport of the measure \x from N to M (the Euler distri- 
bution of mass), and x k are coordinates of the centre of mass in M. For affine 
configurations (I267p implies (|266[) . 

If the internal inertia is isotropic, i.e., (|263l) holds, then (|265[) becomes the 
following geometry-based expression: 

= ^g l3 v l A v^ BV AB = ^G[^ KL n K A n L BV AB 

= ^, J !l , I ff t CM- 1H (268) 

The formula (|262[) works correctly in various problems. Its characteristic 
feature is that it is not invariant under any subgroup of GL(V") x GL(U) larger 
than O (V, j)xO (U, J -1 ) ■ From the purely mathematical point of view this is 
some shortcoming, because on the primary kinematical level, internal degrees of 
freedom are ruled by the total GL(V) x GL(U). Thus, even the very curiosity 
motivated our earlier search of affincly-invariant models of the kinetic energy of 
a single affine body [20l [21] [22] , or in more mathematical terms — the search 
of affinely-invariant Riemannian structures on the affine group. The more so 
that afterwards some possibility of interesting physical applications seemed to 
appear. 

Let us remind the structure of affinely-invariant kinetic energies for a single 
affine body. 

We begin with left-invariant models, i.e., ones invariant under the spatial 
affine group GAf (M). This means they are invariant under (|196p assuming 
that A runs over the total GAf(M), and B = Id^v- If T is to split additively 
into translational and internal parts like in (|262[) . then the only natural and 
at the same time maximally general nonrelativistic models have the following 
form: 

T = T tl + T int = jVabV a V b + ^C B A D c n A B h c D , (269) 

where C B a d c are constants. Obviously, they are symmetric in bi-indices: 

C B A D c = C D C B A- (270) 

Let us observe that such kinetic energies, i.e., metric tensors on the con- 
figuration space Q = M x LI (U, V) are affinely invariant in M, but not in 
the material space N. Namely, (I269P is invariant only under such material 
endomorphisms, i.e., elements of GL(£7), which preserve the material tensors 
ip e U* <g> U* , L £ U <g> U* <g> U <g> U* . The highest possible symmetry in N is 
77-rotational one, i.e., 0(U,rj) C GL(U). It is attained when £ is algebraically 
built of 77 and Idjj, i.e., 

C B a D c = Iv BD Vac + AS b c 6 d a + BS B A 6 D c- (271) 



42 



Then Tj nt in (I269P may be concisely written as 

T int = ~ Tr (fi^fi) + ^Tr (h 2 ) + I f Tr «) * , (272) 



2 V 7 2 V 7 2 
where Q Tv denotes the transpose of the mixed tensor performed in the 77-sense: 



A b=V AC Vbd& D c- (273) 



This is the " usual" transpose when orthonormal coordinates are used in U and 
tjab =* Sab- Obviously, the inertial coefficients /, A, B are constants. 
The metric tensor L underlying f|269[> has the following form: 



mC[tp] tj dx l <g> dx> + -£ A B C Df 1B i<f 1D 3 dtp L A ® <{P b- (274) 



It is evidently curved, (Q, T) is a Riemann space with the non- vanishing cur- 
vature tensor. Compare this with the flat, Euclidean structure of (|264|) . (|265|l , 
where the metric tensor of T was given by 

T = mg lj dx i ® dx j + g lj J AB dip i A ® ip? B . (275) 

Let us observe that the last two terms in (|272[) are affinely invariant both 
in the physical and material spaces. This is impossible for the total kinetic 
energy, when translational motion is taken into account. The combination of 
the mentioned A, £?-terms is not positively definite, for example, the A-term, 
(A/2)Tr (^ 2 )> has the following signature: 

. ( n ( n ^ 1) n(n + 1) \ ,„„„n 
sign A I \ ' -, k 7 > + . (276) 



Roughly speaking, the negative and positive contributions (if A > 0) corre- 
spond respectively to the "compact and non-compact directions" in the linear 
group. The negative contributions to T; nt might seem physically embarrassing, 
nevertheless, they are not only harmless, but even may be desirable in certain 
dynamical models. Incidentally T- mt is positively definite in certain open domain 
of inertial parameters (/, A, B) € M 3 . 

Let us observe that the model (|269j) specialized to the particular case (|27ip , 
1(272]) may be formally "obtained" from the combination of PS3]| , 1(255]) . ((21)51) 
by substituting the constant material metric 77 instead of the ^-dependent 
Green tensor G[ip\. This was just the heuristic motivation for rewriting ([2621) 
in the form (|264[) . (|265p . One obtains then some quadratic function of non- 

holonomic Lie-algebraic velocities (d A , ^b^J with constant coefficients, just as 

(|262p is a constant-coefficients quadratic form of holonomic generalized veloci- 
ties (V, V'Aj- This seemingly formal modification is, as a matter of fact, very 
essential because it leads to affinely- invariant (in M) kinetic energies, i.e., to 
M-affinely metric tensors on the configuration space. 
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Repeating the above way of thinking in terms of Lie-algebraic variables 
(v l ,Q l j) one finds the heuristic way towards models affine-invariant in the 
material space N, i.e., GAff (A^)-invariant ones. The only natural GAff(A)- 
counterpart of (I274p . without interference between translational and internal 
parts has the form: 

T = T tI + T int = yS tf «V + -TVhtfj&u (277) 

with symmetry properties of 1Z like those of C in (|270[) . Dually to (I269p . T tr is 
not invariant under GAff(M), and in a consequence of the structure of the affine 
group, there is no doubly invariant metric on GAff(n,K). However, in analogy 
to (|271[) , there are models invariant under GAff(iV) and simultaneously under 
the Euclidean group ~E(M,g), thus, in particular, under the orthogonal group 
0{V,g). The corresponding formula, dual to ([2"M1) . ([2721 . is 

T int = ^Tr (ft T 9ft) + ^Tr (ft 2 ) + | (Tr ft) 2 , (278) 
where, obviously, 

(n Tg ) 1 3 = g lk g 3 in l k , Tr (ft 2 ) = Tr (ft 2 ) , Tr ft = Tr ft. (279) 
This means that 

TV 'A = Ig jl g ik + AS^S 1 , + BP t S l k . (280) 

Unlike the total kinetic energy, the internal one Ti nt may be simultaneously 
affinely invariant in M and N and is then equal to the sum of the last two 
terms in (|272[) or (|278[) . they are equal to each other. 

It is very instructive to see what would be the most general kinetic en- 
ergy, i.e., metric tensor on the configuration space, invariant under the isometry 
groups in M and N, ordering its terms according to the additional "increasing 
degree" of the affine invariance in space and in the body. 

The corresponding kinetic energy may be written down as follows: 

T = ^( mi G AB +m 2r]AB )d A d B + ^Tv{n 2 ^ +| (Trft) 2 (281) 
+ \ (hG KL G MN + I 2VKLV MN + I 3 G KLV MN + h VKL G MN ) h K M h L N , 
what may be also written in the following alternative form: 

T = l{m 1 g ij +m 2 C ij )v i v j + -Tr(n 2 ) + | (Tr ft) 2 (282) 

+ \ (hg k ig mn + hc kl c mn + hg kl c mn + hc u g mn ) n k m n l n . 

Remark: G KL are not "^-raised" Gab, but their inverses; similarly, C tJ are 
contravariant inverses of C a b, not their contravariant "^-raised" versions. 
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It is very important to say what is the metrical-affine interplay in the ex- 
pression (|28ip or (|282l) . If we put mi, I\, I3, I4, we obtain the model afBncly 
invariant in M but only 77- metrically (and not higher) invariant in N. Models 
g-metrically invariant in M (but not higher invariant there) and affinely invari- 
ant in N are obtained when we choose m%, I2, la, I a- Or, more correctly, there 
is also something between "metrical" and "afBne". And this also has to do 
with the remark a few lines above. Namely, admitting in (|281l) . (|282l) the term 
controlled by I2 , we obtain the model in which the internal kinetic energy is in- 
variant under both 77-rotations and dilatations in the material space, thus, under 
the group of similarities (linear-conformal group R + 0(t/, 77)). Similarly, the I\- 
controlled term is invariant under R + 0(V, g), thus, also under dilatations in the 
physical space. Of course, we mean here the linear-conformal invariance of the 
internal term alone; it does not seem possible to unify the total afhne invariance 
in N with the rotational-dilatational total invariance in M, and conversely. The 
translational term T tr is an obstacle against such a joint invariance. 

Obviously, the metric tensor underlying (I281[) . (12821) 



where 



g = G(tr) tj dx l g> dx 3 + g(mt) A k B idtp k A ® ¥b> (283) 

<7(tr)y = migi j +m 2 Cij, (284) 
0(int) Vi = hg kl G AB + I 2 C klV AB + hg m AB + hC kl G AB 

+ Ay- XA vp- XB k +By- XA W - 1B h (285) 

Let us stress again that £/(int) is autonomous in the sense that its coefficients 
depend only on ip but not on x. Unlike this, £/(tr) is non-autonomous in the 
sense that it depends on internal variables ip (if mi ^ 0). 

It may be also interesting to express (I284p . (|285l) in various non-holonomic 
representations, 

g(ix) = g(tr) AB Q A (286) 
g(tr) AB = miG AB +m 2 r]AB, u A = ip~ 1A t dx* = ip- 1A ,uj\ (287) 
g(int) = g(int) m k n l uj k m ®J n = g(int) B A D cu A B®u c Dl (288) 
£(int)"Y l ; = hgu9 mn + hC kl C mn + hg kl C mn + hC kl g mn 

+ AS m l S n k + BS m k S n h (289) 
g(mt) B A D c = hG AC G BD + I 2 vacV BD + hG Ac v BD + hVA C G BD 

+ AS B c S D A + BS B A S D c, (290) 
u\ - y- lK b d<p a K , i A B = f u #s. (291) 

Roughly speaking, these are metric tensors on the afhne group GAff(n,R) in- 
variant under the left (spatial) and right (material) action of isometries. By 
appropriate specification of constants we obtain metrics on GAff(n, R) with 
larger isometry groups, containing left or right action of afhne transformations. 

This is mathematics, differential geometry, interesting in itself. But there 
are also some physical ideas and hypotheses beyond. Namely, in condensed 
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matter it is quite reasonable to expect that the effective inertial properties of 
molecules are effective, collective resultants of the interaction of any particle 
with its surrounding, instead of being determined by the "true" metric tensor 
of space. There are well-known physical examples in solid state physics, namely, 
the effective mass tensor of the electrons in crystals. 

Our guiding idea is that the search of such effective multiparticle models 
should be based rather on some natural symmetry postulates than on the sys- 
tematic " derivation" from some more " microscopic" structural assumptions and 
calculations, like, e.g., ones on the "atomic" level. In view of the very compli- 
cated system of equations principally describing the problem on the level of 
electrons-nuclei systems, the latter task would be completely hopeless. In such 
situations it happens often that symmetry principles are the only ones which 
enable us to obtain some qualitative and sometimes also quantitative results. 

On the non-dissipative Hamiltonian level, our dynamics will be based mainly 
on the Hamilton functions of the potential shape: 

H = T+V, (292) 

where T is obtained from the kinetic energy models T (metrics on the con- 
figuration space) based on the additive postulate (|258[) . where, for any fixed 
particle K, TQC) has one of the postulated forms (|2"f)2"j) . 1(255)) . (|272)l , ([277]) . 
(|278p . (I28ip . (|282[) . T is obtained from T in a standard way via the Legendre 
transformation. Usually it is convenient to express this transformation in terms 
of non-holonomic velocities like £1 A b, v a , because then T is represented as 
a function of affinc spin in the spatial or co-moving description Y> a b and 
of the linear momentum, also in both possible representations, pi, pa- This is a 
very effective procedure, because equations of motion may be obtained almost 
directly from (|245[) by substituting as phase-space functions F the mentioned 

Hamiltonian generators, i.e., momentum mappings, e.g., (pi, S'j) or (pa, ■ 
The point is that Poisson brackets of these functions are well known, given by 
structure constants of the affine group GAff(n,R). 

For simple models, e.g., (|272[) . (|278[) . the Legendre transformation may be 
explicitly inverted and the first-order Hamilton-type equations (|205|) may be 
easily expressed as second-order differential equations for the time-dependence 
of generalized coordinates x l , f 1 a for the single affine body, and for x l (IC), 
<P % a(IG) for the system of affine bodies. 

Clearly, Poisson brackets for quantities concerning different " particles" , K, ^ 
C, are always vanishing, because they are functions of different generalized co- 
ordinates and momenta. This simplifies remarkably derivation of equations of 
motion. 

Summary 

Let us summarize the ideas and programme outlined here. Discussed was hier- 
archy of dynamical quantities describing systems of affine bodies. In Section 2 
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we reviewed scalar invariants for pairs of affine bodies. The mentioned hierarchy 
concerned the degree of affine invariance. We discussed systems of scalars invari- 
ant under isometries and under the total affine group in physical and material 
spaces. Quite different problem, to some extent experimental one, is whether re- 
ally the affine invariants may be useful in physical models of systems of mutually 
interaction affine bodies. There is also a stronger question: Perhaps in highly 
condensed matter the purely affine scalars might be sufficient for constructing vi- 
able dynamical models of mutual interactions (including both translational and 
internal degrees of freedom) . Is it so or not? There are some physical reasons 
to expect such a possibility. From the purely geometric point of view of rational 
mechanics the problem of such Thales-like physics is interesting in itself. Our 
two-body affine and metrical scalars constructed in Section 2 are thought on as 
argument of binary potentials describing interactions in systems of affine bodies 
("molecules"). Let us repeat: affinely invariant scalars do exist only for systems 
of affine bodies, not for the single affine body. Affinely-invariant kinetic energies 
(metrics on the configuration space) are well defined both for single affine bod- 
ies and for their systems. For the systems we have the additive formulas (|258[) . 
or perhaps the binary ones (|259[) . when something like the gyroscopic coupling 
of internal (angular or affine) velocities occurs. As single kinetic energies, the 
formulas like ([555]) . (12"^ . ([2721 . ([577]) . (1275)) . ([25T1) . ([252]) . etc. may be chose, 
depending on the postulated models, i.e., on the assumed metric tensor of the 
configuration space. Binary expressions T(IC, C) are obtained as so-called po- 
larizations of the one-particle quadratic forms T(/C). The underlying geometry, 
physical motivation and expected usefulness of the particular models of T(/C) 
were discussed in our earlier papers [20l EU [22] . 

The "traditional" model (|262[) is derived for the usual affinely constrained 
system of material points. It is also compatible with the d'Alembert principle 
meant in the sense of the usual Euclidean metric g. In this sense it seems to 
be the most physical model. Nevertheless, there are also some other models, 
like (|269[) - (l28ip . They are geometrically more interesting, because they have 
"large" symmetry groups and the resulting essential nonlinearity in the kinetic, 
inertial part of Lagrangian. And the point is that they seem to offer some 
unexpected dynamical applications. Namely, even on the purely geodetic level, 
without any potential, the doubly affinely invariant models of internal kinetic 
energy, i.e., (|272[) . (|278p without the metrical term, i.e., with I — 0, may encode 
the dynamics of nonlinear elastic vibrations. The general solution contains 
an open family of bounded orbits, even if no potential is assumed. In any 
case, it is so, when the motion is isochoric, i.e., the body is incompressible. 
For the compressible case, some purely dilatational potential may be used to 
stabilized dilatations. And it is interesting that the main characteristics of 
motion my be obtained on the basis of exponential functions of matrices. In 
such bounded solutions, the deformation invariants perform finite vibrations. 
In multibody systems this concerns both the usual deformation invariants and 
invariants of mutual deformation tensors defined in Section 2. The peculiarity 
of multibody affine systems is that one can also use affinely-invariant potential 
energies of binary interactions. This offers a wide class of dynamical models to 
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be postulated as a hypothetical description of interactions in discrete structured 
media. Their continuous limit will lead to a variety of dynamical models of 
structured continua. 

When dealing with objects like molecules, fullerens, graphens and clusters, 
one must use the quantized version of the model. Certain primary ideas in this 
direction were presented in [TSJ [TH] . 
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